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Silicides  produced  from  a codeposited  mixture  of 
gadolinium  and  silicon  have  been  studied  as  an  alternative  to 
the  thin-film  reaction  method.  Results  include  details  of 
the  composition,  crystallographic  structure  and  intrinsic 
properties  of  the  films,  electrical  contact  properties  to 
silicon  and  gallium  arsenide,  and  advancements  in 
quantitative  analysis  of  surface  composition. 

Rutherford  backscattering  spectroscopy  (RBS)  confirmed 
the  reported  temperature  (350-400°C)  of  the  thin-film 
reaction,  and  Auger  electron  spectroscopy  (AES)  depth 
profiles  clarified  the  role  of  oxygen  contamination  in 
inhibiting  the  reaction. 

For  the  codeposited  films  X-ray  diffraction  (XRD)  showed 
progressive  recrystallization  at  temperatures  from  500-800°C, 
with  strong  preferred  orientation  effects  and  evidence  of  the 
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presence  of  the  hexagonal  phase  as  well  as  the  expected 
orthorhombic  phase  in  the  recrystallized  films.  Auger  depth 
profiles  showed  limited  reaction  with  silicon  or  GaAs 
substrates  during  heat  treatment,  with  an  As-rich  compound  at 
the  silicide/GaAs  interface.  Excess  silicon  was  rejected 
from  silicon-rich  films  to  form  the  GdSi1 7 compound. 
Resistivity  of  the  films  was  250^0-cm  after  heat  treatment, 
and  the  Schottky  barrier  height  to  p-Si  was  0.70eV. 

A method  of  deriving  a layered  model  of  near-surface 
composition  is  presented,  using  the  inverse  Laplace  transform 
of  angle-resolved  x-ray  photoelectron  spectroscopy  (ARXPS) 
data.  The  inverse  of  a Laplace  transform  is  not  unique,  and 
thus  constraints  are  required  to  select  a solution  function. 
The  new  approach  relaxes  these  constraints,  allowing  non- 
monotonic profiles  to  be  recovered.  Results  are  presented 
for  an  air-oxidized  Si  substrate  demonstrating  a layered 
surface  oxide,  with  silicon-rich  oxides  at  the  Si02/Si 
interface.  Ar-sputtered  GaAs  showed  a surface  enrichment  of 
As  above  a Ga-rich  layer,  consistent  with  preferential 
sputtering  by  Gibbsian  surface  segregation  of  As. 

Empirically  derived  AES  correction  factors  are  presented 
that  permit  accurate  quantification  of  rare  earth  silicide 
films,  correcting  for  effects  ranging  from  surface 
segregation  to  changes  in  the  effective  attenuation  length  of 
the  Auger  electrons  and  backscatter  effects. 
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CHAPTER  ONE 
INTRODUCTION 

Purpose 

Large  scale  implementation  of  photovoltaic  power 
generation  is  dependent  to  a great  extent  on  advances  in 
materials  and  processing  to  achieve  greater  solar  cell 
efficiency.  The  incentive  is  particularly  strong  in  solar 
concentrator  applications,  where  the  marginal  economic 
return  is  multiplicative  compared  to  the  cost  of  increased 
cell  capacity. 

The  purpose  of  this  study  has  been  to  evaluate  the 
potential  of  rare  earth  silicides  for  improving  solar  cell 
performance.  The  role  of  the  silicides  would  be  to  serve  as 
a contact  metallization  layer,  by  which  electrical  current 
is  passed  into  or  out  of  the  semiconductor  solar  cell 
material.  Improvements  in  contact  metallization  can  enhance 
concentration  cell  performance  by  reducing  direct  losses  due 
to  contact  resistance  and  surface  recombination  of  carriers. 
In  addition,  the  service  lifetime  of  the  cells  may  be 
extended  by  enhanced  metallurgical  stability  of  the  metal- 
semiconductor  contact. 

This  means  that  an  understanding  of  the  electrical 
transport  properties  of  the  silicide  and  of  the  silicide- 
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semiconductor  interface  are  of  central  importance  to  this 
study.  To  some  extent  these  properties  may  be  described  by 
measurements  of  such  parameters  as  the  resistivity  or  thin 
film  sheet  resistance  of  the  silicide,  or  the  Schottky 
barrier  height  or  specific  contact  resistance  of  contacts. 

However,  a full  and  complete  characterization  of  the 
samples  is  as  necessary  to  the  study  as  the  measurement  of 
such  device  performance  parameters.  The  actual  sample 
morphology,  the  phases  present,  condition  of  surfaces  and 
interfaces,  and  a host  of  other  attributes  may  be  expected 
to  vary  considerably  with  variations  in  preparation  and 
treatment  and  can  affect  the  response  of  samples  to 
measurement. 

Analysis  of  the  deposited  films  has  required  the 
application  of  numerous  techniques,  and  reconciliation  of 
the  results  of  the  various  techniques  to  give  a full  picture 
of  the  material  under  study.  Fuller  understanding  of  the 
data  led  to  some  of  the  significant  results  here  presented, 
including  the  development  of  novel  techniques  in  data 
reduction  and  corrections  to  standard  analyses. 

Overview 

The  effort  to  place  this  work  properly  in  its  context 
begins  with  a review  of  the  literature.  Chapter  Two 
proceeds  from  a simple  model  of  the  metal-semiconductor 
junction  to  metallization  requirements  of  semiconductor 
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devices  in  general,  with  models  and  analyses  of  several 
performance  parameters,  as  an  outline  of  the  current 
technology.  The  chapter  continues  with  a review  of  the 
basic  properties  of  thin  films  of  other  metal  silicides  as  a 
background  for  an  introduction  to  previous  work  with  rare 
earth  silicides  in  the  form  of  reacted  thin  films  on 
silicon. 

Chapter  Three  describes  experimental  work  including 
fabrication  of  samples  and  fundamental  thin  film  analysis  to 
identify  the  disilicide  phase.  Early  work  involving  the 
deposition  of  rare  earth  metals  for  reaction  with  silicon 
substrates  provided  samples  for  comparison  with  codeposited 
films,  permitting  the  development  of  the  ion  beam  sputter 
codeposition  technique  for  gadolinium  silicides. 

Chapters  Four  and  Five  present  examples  of  two  areas  in 
which  a more  in-depth  approach  to  data  analysis  has  been 
developed  and  demonstrated  as  necessary  to  enhance  the 
accuracy  of  two  commonly  used  analytical  techniques  for 
quantitative  analysis  of  surfaces.  Chemical  composition  as 
measured  by  electron  spectroscopic  analysis  of  a clean 
surface  will  often  deviate  from  the  bulk  composition  of  a 
substrate  or  thick  film.  In  Chapter  Four  it  is  shown  that 
changes  in  near-surface  composition  may  have  a significant 
effect  on  the  angular  distribution  of  photoelectrons 
escaping  the  sample,  with  resultant  nontrivial  variations  in 
the  calculated  relative  abundance  of  the  sample's  elemental 
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constituents.  It  is  further  shown  that  from  such 
measurements,  taken  at  a series  of  angles  of  observation 
with  respect  to  the  sample  surface,  it  is  possible  to  derive 
layered  models  of  the  composition  changes  in  a shallow 
region  near  the  sample  surface.  In  Chapter  Five,  correction 
factors  are  derived  to  permit  accurate  guantif ication  of 
rare  earth  silicide  films  using  Auger  electron  spectroscopy 
(AES) . These  empirically  derived  factors  correct  for 
effects  ranging  from  surface  segregation  to  changes  in  the 
effective  attenuation  length  of  the  Auger  electrons  and 
backscatter  effects. 

In  Chapter  Six  are  the  detailed  results  of  analyses. 
Electrical  contacts  to  both  silicon  and  gallium  arsenide 
have  been  studied,  as  well  as  the  intrinsic  physical, 
electrical  and  optical  properties  of  the  silicides.  Strong 
crystallographic  orientation  effects  were  observed  in  the 
recrystallized  films  and  are  believed  to  play  an  important 
role  in  contact  properties. 

Chapter  Seven  presents  a summary  of  the  work  and 
concluding  remarks,  with  observations  relevant  to 
applications  of  these  films  and  suggestions  for  future 
studies. 

The  balance  of  this  chapter  introduces  the  metal 
silicides  as  semiconductor  contacts,  and  the  rare  earth 
elements  and  their  silicides  as  interesting  candidates  for 
forming  high-temperature  stable  Ohmic  contacts. 
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Silicides  for  Semiconductor  Metallization 

Metal  silicides  have  recently  emerged  as  important 
contact  materials  for  semiconductor  applications,  mainly  as 
gates  and  interconnects  in  integrated  circuit  technology. 
Most  of  the  progress  has  been  realized  using  silicides  of 
the  transition  metals,  that  is  those  elements  in  which 
valence  electrons  are  being  added  to  the  3d,  4d  and  5d 
orbitals  across  the  periodic  chart.  These  silicides  can  be 
formed  by  reaction  of  a deposited  metal  layer  with  a silicon 
substrate  at  elevated  temperatures.  The  silicide  films 
formed  in  this  reaction  have  advantages  of  good  stability  at 
high  temperatures  and  stable  interfaces  with  silicon. 

Silicide  films  can  also  be  formed  by  codeposition  with 
annealing  to  form  the  silicide.  Codeposited  silicide  films 
have  found  application  in  silicon  technology  and  in  GaAs 
devices. 

The  properties  of  transition  metal  silicides  are 
summarized  in  Chapter  Two.  By  virtue  of  the  kinetics  of 
their  formation  and  properties,  these  metal  silicides 
naturally  fall  into  two  classes,  the  refractory  metal 
silicides  (metals  from  groups  IVB  to  VIB  in  the  periodic 
chart)  and  the  silicides  of  the  noble  and  near-noble  metals 
(group  VIII) . 

The  rare  earth  silicides  form  in  a reaction  that  is 
highly  exothermic  and  in  many  ways  unigue  among  the  metal- 
silicon  systems.  They  constitute  an  additional  extensive 
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class  of  silicides  with  different  reaction  mechanisms, 
crystal  structure  and  film  morphology.  However,  their 
properties  have  been  far  less  investigated. 

The  Rare  Earths 

The  rare  earth  metals  (i.e.,  the  elements  in  the 
lanthanide  seguence,  from  La  to  Lu)  and  their  compounds 
often  exhibit  unusual  properties1,2  due  to  the  uncommon 
electronic  structure  of  these  elements,  with  electrons  being 
added  to  an  inner  4f  shell  that  is  lower  in  energy  than  the 
bonding  orbitals,  6s  and  5d.  These  inner  shell  4f  electrons 
are  well  shielded  and  tend  to  be  localized  to  their 
respective  atoms.  Rather  than  being  continuously  modified 
by  the  addition  of  valence-level  electrons  available  for 
bonding,  lanthanide  chemical  properties  are  to  a large 
extent  guite  similar  across  the  seguence.  Chemical  bonding 
involves  only  the  6s  and  5d  electrons  common  to  all  the 
elements,  while  the  distinguishing  electron  from  element  to 
element  is  added  to  the  more-tightly  bound  4f  shell. 

This  unusual  electronic  structure  is  responsible  for 
unusual  behavior  such  as  the  Lanthanide  contraction,  the 
decrease  in  atomic  radius  with  the  increasing  effective 
nuclear  charge  across  the  series.  These  elements  are 
comparable  to  magnesium  in  electropositive  character  and 
ease  of  oxidation.  Removal  of  the  valence  electrons  in 
compound  formation  produces  a +3  oxidation  state  that  is  the 
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normal  state  for  all  of  the  lathanides,  with  the  outermost 
remaining  electrons  being  the  distinguishing  4f  shell. 

Excitation  of  these  4f  electrons  requires  considerable 
energy,  and  optical  absorption  in  the  visible  range  for  rare 
earth  compounds  is  unusual  in  that  it  tends  to  occur  in 
sharp  peaks  similar  to  gas  phase  absorption  rather  than  the 
broadened  bands  typical  for  condensed  states.  Mechanisms  of 
electron  transport  by  hopping  between  localized  states  are 
significant  to  the  total  conductivity  of  these  compounds  as 
well . 3 

Bonding  is  generally  trivalent,  involving  the  6s  and  5d 
orbitals  that  they  all  possess  at  the  valence  level.  In 
addition  to  the  lanthanides,  yttrium  and  scandium  are  often 
included  under  the  heading  of  "rare  earths"  because  of  their 
similarity  to  lanthanum  as  precursors  in  the  IIIA  column  of 
the  periodic  chart. 

High  Temperature  Ohmic  Contacts 

It  is  often  desirable  to  inject  current  into  or  out  of 
a semiconductor  with  a low  power  loss.  Such  a contact 
should  ideally  present  a negligible  barrier  to  current  and 
thus  should  obey  Ohm's  law  (the  current  should  increase 
linearly  with  increasing  voltage) . The  constant  of  this 
proportionality  is  a resistance  attributed  to  the  contact 
interface. 
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Metal  silicides  are  generally  known  to  be  stable  at 
high  temperatures  and  can  make  ohmic  contacts  to 
semiconductors.4'5  Previous  work  with  reacted  thin  films  of 
rare  earth  metals  on  silicon  substrates  has  demonstrated 
that  rare  earth  metal  films  will  react  with  silicon  at  low 
temperatures  (300-400C0)  to  form  stable  disilicides  with 
extremely  low  Schottky  barrier  heights  to  n-Si.6,7 

Ohmic  contacts  to  n-type  silicon  are  predictable  from 
these  low  Schottky  barrier  heights,  while  contact  properties 
to  GaAs  using  these  silicides  have  not  been  reported.  Rare 
earth  metals  like  other  metals  have  higher  Schottky  barrier 
heights  to  n-GaAs,  and  silicide  contacts  may  be  expected  to 
depend  upon  enhanced  tunnelling  or  field  emission  to  produce 
ohmic  behavior. 

Successful  ohmic  contacts  to  n-GaAs  have  rarely 
involved  simple  tunneling  behavior  through  planar 
interfaces.  Alloyed  contacts,  in  which  reaction  with  and 
consumption  of  portions  of  the  GaAs  substrate  are  believed 
to  create  a highly-doped  surface  layer,  produce  a nonplanar 
and  inhomogeneous  interface  with  alloy  penetrations  which 
increases  the  contribution  of  field  emission  to  the  current 
flow,  as  well  as  increases  the  effective  area  of  the 
contact.8 

Such  an  approach  is  limited  at  higher  temperatures  if 
continued  diffusion  and  reaction  lead  to  degraded  contact 
properties  or  device  failure.  Rare  earth  silicides  have  the 
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potential  for  forming  more  abrupt  and  thermally  stable 
interfaces  to  GaAs , while  their  ability  to  alloy  and 
transfer  dopants  to  the  substrate  may  be  limited.  Thus  they 
must  depend  heavily  upon  low  barrier  height  and  planar 
tunnelling  for  forming  ohmic  contacts  to  GaAs. 

Gadolinium  as  a Test  Case 

The  Lanthanide  sequence  has  unique  properties  that  tend 
to  follow  orderly  trends  across  the  periodic  chart,1  and  if 
a property  of  one  rare  earth  falls  in  an  undesired  range, 
one  has  reasonable  hope  that  further  up  or  down  the  series 
one  may  find  an  element  that  works.  Gadolinium  presents  a 
good  first  choice  for  study,  falling  in  the  middle  of  the 
sequence.  It  has  a 4f  electron  shell  that  is  half-full  with 
7 electrons,  which  is  a configuration  of  enhanced  stability 
that  affords  some  unique  properties  among  the  lanthanides. 

A search  of  the  literature  yields  somewhat  more  information 
on  the  Gd/Si  system  than  the  other  rare  earth  silicides. 

The  choice  of  Gd  carries  with  it  some  advantages  and 
some  disadvantages.  The  lighter  lanthanides  (Z<64)  are  less 
attractive  choices  due  to  their  higher  reactivity  with 
oxygen.  However,  some  of  the  heavier  lanthanides  may  offer 
greater  potentials  in  that  the  hexagonal  phase  becomes  more 
stable.  There  is  evidence  that  this  hexagonal  phase  may 
offer  enhanced  conductivity  based  on  the  results  of  Thompson 
et  al.9  This  is  a reasonable  supposition  in  view  of  the 
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metal/silicon  layered  structure  in  the  hexagonal  phase.  As 
the  metal  atoms  decrease  in  size,  the  possibility  of 
epitaxial  growth  of  the  hexagonal  phase  is  also  enhanced  due 
to  improved  lattice  matching  with  the  Si (111)  face.  This 
match  is  optimized  for  ErSi2.10 

However,  gadolinium  falls  on  the  borderline  as  far  as 
stability  of  the  orthorhombic/hexagonal  phases,  and  has  a 
demonstrated  potential  to  grow  epitaxially  on  either  the 
(100)  face11,12  or  the  (111)  face13  of  silicon.  Among  the 
rare  earth  silicides,  it  has  the  lowest  barrier  height  to 
n-type  silicon.6,7  Gadolinium  was  therefore  chosen  as  a 
particularly  interesting  element  to  study  as  a 
representative  of  the  rare  earth  sequence. 


CHAPTER  TWO 
REVIEW  OF  LITERATURE 

Metallization  Requirements 

The  role  of  metallizations  can  be  seen  in  four  basic 
applications.  At  the  metal  level  in  integrated  circuit 
devices  there  are  two  functions,  as  gate  or  as  interconnect 
metallizations.  The  MOSFET  (metal  oxide  semiconductor  field 
effect  transistor)  is  controlled  by  a gate  voltage  which  is 
applied  to  a metal  electrode  that  is  isolated  from  the 
substrate  by  a dielectric  layer.  This  electrode  and  other 
device  terminals  are  interconnected  by  a network  of  metal 
runners,  called  the  device  interconnect  metallization. 
Current  state-of-the-art  devices  use  polysilicon  (with  a 
sheet  resistance  near  30-60  ohm/sq.)  in  both  roles.14 

In  addition,  current  must  enter  or  leave  the 
semiconductor  portions  of  a device  through  contacts,  or 
metal-semiconductor  interfaces,  which  can  function  as  either 
a barrier  (current-rectifying)  contact  or  as  an  ohmic . low 
resistance  (current-passing)  contact. 

The  interest  in  silicides  is  based  on  the  need  to 
develop  new  metallization  schemes  for  the  microelectronics 
industry.  Gates  and  interconnects  are  subject  to 
increasingly  stringent  requirements  as  VLSI  (very  large 
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scale  integrated  circuit)  developments  lead  to  devices  that 
are  smaller  in  size,  faster  in  response  and  consuming  less 
power.  Intrinsic  RC  delay  due  to  high  sheet  resistance  has 
placed  a severe  limit  on  the  usefulness  of  polysilicon  films 
as  the  conducting  metallization. 15  As  device  size  is 
reduced,  fringe  effects  on  gate  capacitance  begin  to 
dominate  and  it  becomes  clear  that  the  only  potential  for 
reducing  the  RC  time  constant  is  to  reduce  the  sheet 
resistance  term.14 

On  the  other  hand,  most  metals  are  subject  to 
limitations  in  their  stability  under  device  processing  or  in 
the  operating  ambient.  Whether  the  problem  is  low 
temperature  reactions  with  the  substrate  or  dielectric 
layer,  low  temperature  melting  or  eutectic  formation,  high 
oxidation  rate  or  susceptibility  to  electromigration,  the 
list  of  choices  among  the  low-sheet-resistance  metals 
dwindles  rapidly.14  Thus  a prime  motivation  is  to  provide 
alternative  low-sheet-resistance  conductors  to  enable  the 
smaller  device  structures  to  reach  their  potential  for 
faster  switching  speeds. 

Contact  applications  are  also  affected  by  the  size 
reduction  in  VLSI  developments.  Reduced  junction  depth 
means  increased  problems  of  metal  penetration  into  silicon, 
especially  for  aluminum.14  Silicides  formed  by  reaction  may 
have  a similar  disadvantage  in  that  silicon  is  consumed  in 
the  reaction  and  shallow  junctions  may  therefore  be 
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penetrated.  But  silicides  that  are  codeposited  as  a mixture 
can  minimize  the  substrate  reaction,  and  are  expected  to  be 
more  stable  over  shallow  junctions . 16,17,18 

Metal-Semiconductor  Contacts 
Schottkv  Barriers 

The  simplest  model  of  a metal-semiconductor  interfacial 
barrier  to  current  is  the  Schottky  barrier  model.  In  this 
model  the  Fermi  energy  of  the  metal  and  the  semiconductor 
are  equalized  by  a transfer  of  charge  at  the  interface, 
resulting  in  a dipole  electric  field  at  the  interface.19 
The  metal  carries  a charge  at  its  surface,  while  the 
semiconductor  distributes  an  equivalent  but  opposite  charge 
across  a finite  layer  known  as  a space  charge  region, 
equivalent  to  that  formed  at  a semiconductor  n-p  junction. 
Since  the  majority  carriers  are  opposed  to  the  charge 
developed,  they  are  depleted  in  number  within  the  space- 
charge  region,  which  is  therefore  also  referred  to  as  the 
"depleted"  region. 

The  charge  imbalance  at  the  interface  acts  as  a 
potential  barrier  to  the  flow  of  carriers  across  the 
interface,  and  a rectifying  contact  results.  Charge 
transport  under  forward  bias  is  governed  by  the  process  of 
thermionic  emission,  with  only  those  carriers  having 
adequate  thermal  energy  to  cross  the  barrier  available  to 
carry  the  current.  The  current  density  across  a Schottky 
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barrier  of  height  <pB  with  applied  voltage  V is  governed  by 
the  equation 

J = A**T2  [exp  ( —j^r ) 1 [exp(~^'1)]  (2-1) 


where  A**  is  Richardson's  constant,  k is  Boltzmann's 
constant,  T is  the  temperature  (Kelvin) , and  n is  the 
"ideality"  factor,  which  is  unity  for  ideal  Schottky 
contacts. 

According  to  this  equation  a plot  of  the  logarithm  of 
the  current  versus  applied  voltage  should  be  linear,  with 
the  barrier  height  derivable  from  the  zero-voltage 
intercept.  This  approach  is  commonly  used  to  measure 
Schottky  barrier  heights,  assuming  an  abrupt,  planar 
interface  between  a semiconductor  substrate  and  a uniform 
metal  layer.  However,  even  with  elemental  metal  films  on 
elemental  substrates  this  is  not  always  the  case  and  is 
increasingly  unlikely  with  compound  films  and/or  compound 
substrates.20 

Schottky  barrier  height  and  work  function 

In  the  pure  Schottky  model,  the  Schottky  barrier  height 
is  defined  as  the  difference  between  the  electron  affinity 
of  the  semiconductor  and  the  work  function  of  the  metal. 

The  work  function  of  a material  is  difficult  to  measure 
accurately,  but  it  has  been  pointed  out  that  the  work 
function  is  related  to  the  electronegativity  of  a metal,21 
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which  is  a measure  of  its  tendency  to  acquire  an  additional 
electron.  Plots  of  barrier  height  against  either  quantity 
show  reasonable  agreement  with  the  theory.  For  a given 
elemental  or  compound  semiconductor,  barrier  height  is 
dependent  upon  the  ionicity  of  the  bonds  in  the 
semiconductor . 22 
The  Bardeen  model 

Real  contacts  rarely  exhibit  ideal  Schottky  behavior. 
Surface  electron  states  due  to  defects  (such  as  unsatisfied 
bonds  in  the  terminated  semiconductor  lattice  or 
contamination  of  the  interface)  can  produce  departures  from 
the  expected  or  "ideal"  barrier  height.23  These  surface 
states,  occurring  at  forbidden  energies  in  the  semiconductor 
band  gap,  are  capable  of  compensating  the  charge  of  the 
metal,  which  in  turn  pins  the  Fermi  level  of  the 
semiconductor  surface  at  the  level  of  the  filled  traps. 

This  results  in  a barrier  height  that  is  independent  of  the 
metal  work  function.  While  for  silicon  either  model  has 
application,  for  the  case  of  compound  semiconductors  such  as 
GaAs , the  Bardeen  model  is  necessary  to  explain  this 
observed  independence  of  barrier  height  from  the  metal  work 
function. 

Schottky  barriers  of  metal  silicides 

For  the  reacted  silicides,  the  scavenging  effect  of  the 
silicide  formation  leads  to  an  interface  that  is  relatively 
free  of  imperfections  and  contamination.  For  this  reason, 
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reliable  and  reproducible  Schottky  barrier  contacts  are 
readily  achieved  using  silicide  reactions.24  Devices  with 
invariant  contact  properties  and  a range  of  barrier  heights 
are  produced  in  this  fashion. 

Barrier  heights  for  the  silicides  are  not  far  in  value 
from  those  of  the  metals  themselves  and  have  been  plotted 
relative  to  the  metal  work  function  with  reasonable 
agreement  to  linearity,  although  the  silicon-rich  silicides 
tend  to  have  lower  barrier  heights  than  other  silicides  and 
less  dependence  on  the  work  function  of  the  metal.25  Plots 
against  a weighted  average  work  function  using  an  arbitrary 
hypothetical  interfacial  composition  do  not  yield  improved 
results.26 

The  relation  of  barrier  height  to  work  function  is  not 
totally  satisfactory,  and  other  correlations  are 
observable.5  Many  transition  metals  sharing  a high  affinity 
for  oxygen  (groups  IVB-VIB)  exhibit  a barrier  height  near 
half  the  silicon  band  gap.  Metals  in  group  VIII  have  little 
affinity  for  oxygen  and  barrier  heights  that  are  high  and 
increasing  with  the  number  of  d-electrons. 

Other  approaches  taken  involve  thermodynamic 
considerations  that  can  be  related  to  electronegativity  and 
the  work  function,  such  as  the  heat  of  formation  of  the 
silicide,27  or  eutectic  temperatures  from  the  binary  phase 
diagrams.28  The  first  approach  predicts  no  barrier  higher 
than  0.83 eV , which  fails  for  PtSi  and  IrSi.24  The  latter 
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approach  involves  choosing  the  temperature  of  that  eutectic 
composition  (in  the  silicon-metal  binary  phase  diagram) 
which  is  closest  to  the  side  of  the  dominant  diffusing 
species  in  the  silicide-forming  reaction  and  is  reasonably 
successful.  It  is  postulated  that  the  properties  of  the 
reactive  interface  are  dominated  by  an  interfacial  layer  of 
a eutectic  composition,  which  is  energetically  favored  to 
form  when  the  appropriate  composition  is  reached  at  the 
interface  due  to  the  interdiffusion  of  the  two  species.  It 
has  been  suggested  that  a low  Schottky  barrier  height  is 
related  to  the  eutectic  temperature  of  the  silicide,  with  a 
decreasing  barrier  height  associated  with  increasing 
eutectic  temperature. 

Ohmic  Contacts 

The  simplest  model  for  the  establishment  of  an  ohmic 
contact  to  a semiconductor  is  the  Schottky-limited  model, 
where  the  contact  is  a simple  junction  between  the 
semiconductor  and  a metal  that  has  a work  function  closely 
matched  to  the  semiconductor  conduction  band-edge  (for 
n-type) , or  the  valence  band-edge  (for  p-type) . A work 
function  slightly  above  (or  below)  the  band  edges  of  the 
semiconductor  causes  an  accumulation  rather  than  a depletion 
of  carriers,  which  presents  no  barrier  to  current  flow  for 
either  the  n-type  (or  the  p-type)  semiconductor.29  A metal 
having  a work  function  within  the  band  gap,  but  near  a band 
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edge,  has  a very  low  Schottky  barrier  height  to  that 
particular  band,  with  the  contact  closely  approaching  ohmic 
behavior. 

In  actuality  it  is  found  that  metals  of  suitable  work 
function  are  seldom  available  or  practical  in  application, 
even  if  Fermi  level  pinning  were  not  found  to  play  a 
dominant  part  in  barrier  creation.  However,  ohmic  behavior 
can  still  be  approached  if  a high  carrier  density  is  present 
at  the  surface  of  the  semiconductor,  thus  creating  a very 
narrow  space-charge  region  which  can  be  penetrated  by 
tunneling  of  the  carriers.30 

Thus  formation  of  practical  ohmic  contacts  of  low 
contact  resistance  is  generally  achieved  by  forming  a highly 
doped  layer  in  the  semiconductor  in  contact  with  the 
metal.31  Higher  doping  levels  in  the  semiconductor  cause  a 
narrowing  and  lowering  of  the  potential  barrier.  The  space 
charge  region  is  made  narrower  because  the  increased  dopant 
concentration  allows  charge  balance  to  be  achieved  over  a 
shorter  distance  (expressed  as  shortening  of  the  Debye 
length,  the  characteristic  distance  over  which  a charge 
imbalance  is  dissipated  in  the  semiconductor) . In  addition, 
the  barrier  height  is  lowered  due  to  image  force  effects 
(charge  redistribution  when  a carrier  approaches  the 
barrier) . The  narrower  and  lower  the  barrier,  the  greater 
the  probability  of  quantum  mechanical  tunneling,  and 
increased  tunneling  of  carriers  is  often  used  to  produce 
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contacts  that  are  effectively  ohmic,  with  tunneling  as  the 
dominant  conduction  mechanism. 

This  is  practically  realized  by  a number  of  approaches. 
Some  depend  upon  the  formation  of  a low-melting  eutectic 
interfacial  phase  in  a metallurgical  reaction  that  consumes 
portions  of  the  semiconductor.  In  some  cases  a highly  doped 
surface  layer  of  the  semiconductor  recrystallizes  from  the 
molten  eutectic  (e.g.,  Al/Si) .32  Others  depend  on  the 
indiffusion  of  a dopant  from  the  metal  layer  during  the 
reaction.  In  the  case  of  n-GaAs,  ohmic  contacts  have  been 
achieved  using  a Au-Ge  or  Ni-Au-Ge  alloy.  It  is 
hypothesized33  that  a liquid  phase  eutectic  is  formed  at  the 
interface  during  heat  treatment,  with  Ga  outdif fusion  into 
the  liquid  permitting  Ge  to  penetrate  the  GaAs  lattice  and 
occupy  the  empty  Ga  sites.  This  creates  an  n+  surface  layer 
and  a sharply  defined  space-charge  region,  with  tunneling 
through  the  barrier  as  the  major  mechanism  of  charge 
transport  at  the  interface. 

Devices  using  contacts  that  are  based  on  reactive 
metallurgy  often  fail  due  to  excessive  reaction  between  the 
semiconductor  and  the  metal  layer.34,35  Whether  it  takes  the 
form  of  "spiking"  (penetrations  of  reacting  metal  into  the 
semiconductor,  sometimes  to  the  cell  junction  causing  a 
junction  short) , or  recombinative  poisoning  by  electrically 
active  metal  atoms  diffusing  into  the  semiconductor,  these 
contacts  are  subject  to  failure  in  processing  and  in  the 
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field.  The  very  reactivity  that  activates  these  contacts  is 
capable  of  corrupting  them.  The  dependence  on  these 
reactions  can  be  avoided  by  preparing  a substrate  with  a 
highly  doped  ion-implanted  or  grown  surface  layer,  and  using 
a metal  that  is  unreactive  with  the  substrate. 

Special  application — solar  cells 

In  developing  a solar  cell  metallization  scheme,  there 
are  a handful  of  criteria  that  can  be  decisive  in  satisfying 
the  question  of  its  usefulness:  1.  acceptable  cost,  2.  low 
contact  resistance,  3.  low  sheet  resistance,  4.  stability  in 
processing  (including  heat  treatments),  and  5.  durability  in 
service. 

For  concentrator  applications,  a higher  cost  may  be 
acceptable,  but  a lower  contact  resistance  and  a lower  sheet 
resistance  of  the  metal  are  desired  because  of  the  higher 
current  densities  that  must  be  handled  by  the  cell 
metallization.36  Thus,  multilayer  structures  are  commonly 
used,  with  a contact  layer  providing  a good  ohmic  contact  to 
the  semiconductor  as  well  as  a diffusion  barrier  that  should 
provide  good  adhesion  for  a thick  conducting  overlayer. 

Metal  Silicides  for  ohmic  contacts 

With  a stable  interface  to  silicon,  and  assurance  of 
reduced  junction  penetration  by  using  codeposited  silicides 
or  silicon-metal  mixtures,  silicides  have  the  potential  for 
forming  contacts  of  low  resistance,  particularly  in  the  case 
of  the  near-noble  metal  silicides.37,38 


The  barrier  height 
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can  be  reduced  dramatically  by  controlling  the  surface 
doping  of  the  semiconductor,19  and  ohmic  contacts  should  be 
readily  achievable  for  most  of  these  silicides,  especially 
the  noble  metal  silicides  on  p-silicon. 

Metal  silicides  exhibit  several  properties  that 
recommend  them  as  candidates  for  forming  stable  ohmic 
contacts;  they  have  good  high  temperature  stability  and  can 
form  stable  interfaces  of  low  contact  resistance  to  silicon. 
They  have  further  potential  application  as  contacts  to 
gallium  arsenide  if  their  reactivity  with  the  gallium 
arsenide  is  limited. 

Contact  Resistance 

Once  an  ohmic  contact  is  achieved,  a criterion  of 
quality  is  a low  resistance  to  the  passage  of  current. 
Specific  contact  resistance  (ohms-cm2)  is  defined  as  the 
resistance  of  a unit  area  of  the  metal-semiconductor 
interfacial  layer.39  Contact  resistance  is  a strong 
function  of  the  doping  level  of  the  semiconductor,  which  is 
reasonable  in  that  the  Schottky  barrier  height  and  Debye 
length  are  sensitive  to  the  carrier  concentration. 19  Thus 
there  is  great  potential  in  processing  for  optimizing  the 
contact  resistance  of  tunneling  ohmic  contacts. 

The  transmission  line  model  (TLM)  of  the  metal- 
semiconductor  contact  introduced  by  Shockley  and  refined  by 
Berger40  is  used  in  a common  method  of  measuring  specific 
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contact  resistance.  Measurements  on  a test  structure 
consistent  with  the  model  permit  the  evaluation  of  a 
transfer  length,  which  in  effect  expresses  the  contribution 
of  the  contact  resistance  to  the  total  resistance  of  the 
test  structure  as  equivalent  to  an  additional  length  of 
substrate. 

Contact  resistance  may  also  be  given  in  units  of  line 
resistance  (ohm-mm) . This  "front"  resistance  gives  the 
interfacial  resistance  in  terms  of  the  width  of  the  contact 
and  is  equivalent  to  the  product  of  a transfer  length  and  a 
sheet  resistance  of  the  substrate.  It  is  evaluated  directly 
from  the  measurements  and  thus  incorporates  any  changes  in 
the  substrate  resistivity  due  to  formation  of  the  contact. 
Contact  resistance  measurement 

The  transmission  line  model  (TLM)  allows  the 
calculation  of  specific  contact  resistance  based  upon 
measurements  of  voltage  and  current  across  a suitable  test 
structure.40  The  structure  generally  used  consists  of 
rectangular  contact  bars  of  width  W,  with  a length  L and  a 
spacing  S along  the  current-carrying  axis.  Conduction  is 
restricted  to  a thin  diffused  or  doped  surface  layer  of 
thickness  t (t=0. 3-1.0  microns),  by  creating  a junction 
potential  barrier  that  prevents  carriers  in  the  bulk  of  the 
substrate  from  participating  in  the  charge  transport.  The 
model  requires  that  S»t,  W»L,  and  that  current  spreading 
beyond  the  pattern  be  negligible. 
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Correct  measurements  to  accurately  separate  the  contact 
resistance  from  the  other  resistances  in  the  TLM  structure 
are  not  always  straightforward.  Accurate  measurement  of 
very  low-resistance  contacts  is  difficult.  Probe  placement 
may  have  an  undesirable  effect,  particularly  for  the  case  of 
a more  highly  resistive  film,  and  the  "front"  and  "end" 
contact  resistances  may  be  confused.  Measurements 
necessarily  include  both  forward-  and  reverse-biased 
contacts,  which  for  marginal  contacts  can  have  the  effect  of 
exaggerating  the  forward  resistance  of  interest  for  solar 
cell  applications.  Also,  the  sheet  resistance  of  the 
substrate  may  be  altered  beneath  the  contacts  due  to 
reactions  in  forming  the  contact  (e.g.,  scavenging  some 
thickness  of  the  semiconductor,  doping  by  diffusion, 
creating  Auger  recombination  centers) , the  effects  of  which 
are  not  accounted  for  in  this  model. 

Other  test  structures 

Meier  and  Schroder41  presented  another  structure 
consisting  of  an  array  of  narrow  parallel  bars  (n=40) , in 
which  the  effects  of  contact  resistance  are  accumulated  over 
the  pattern,  which  reduces  the  effects  of  current  and 
voltage  probes.  In  addition,  this  model  allows  for  the 
sheet  resistance  of  the  substrate  under  the  contact  bars  to 
be  different  from  the  sheet  resistance  between  the  bars. 

A four-terminal  Kelvin  resistance  test  structure39,42 
allows  the  direct  measurement  of  interfacial  contact 
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resistance,  obviating  the  need  to  evaluate  a transfer  length 
which  is  small  compared  to  the  other  terms  involved. 
Limitations 

Modeling  the  current  flow  in  various  test  structures 
shows  that  contact  resistance  measurements  may  vary 
considerably  due  to  a variety  of  factors,  some  geometric  and 
some  related  to  sample  properties;43  attaching  a clear 
meaning  to  a particular  measurement  is  not  at  all 
straightforward . 

The  test  structure  should  resemble  the  device  geometry 
in  order  for  the  measurement  to  be  an  appropriate  assessment 
of  contact  resistance  for  that  application.  The  Kelvin  test 
structure  referred  to  above  is  capable  of  determining  very 
low  contact  resistances,  but  it  also  measures  a current 
passing  through  a plane  of  contact  and  thus  reduces  edge 
effects.  However,  for  the  application  at  hand  as  a top- 
contact  for  solar  cells,  current-crowding44  at  the  edge  of 
the  contact  due  to  the  predominately  horizontal  flow  of  the 
current  is  expected  to  dominate  the  device.  In  addition, 
since  losses  due  to  contact  resistance  are  expected  to  be  on 
the  order  of  0.5%  of  the  power  loss,  practical  contact 
resistances  need  not  be  as  low  as  those  that  necessitate  the 
Kelvin  measurement. 

Measurements  using  the  TLM  model  depend  upon  the 
conductance  of  the  metal  film  being  very  much  greater  than 
that  of  the  substrate,  and  furthermore  that  the  current 


25 


carried  in  substrate  be  confined  to  a thin  surface  layer, 
modelled  as  a transmission  line.  For  more  resistive  metals, 
the  first  requirement  (i.e.,  that  the  conductance  be  higher 
than  that  of  the  substrate)  would  call  for  a thinner  or  more 
lightly  doped  conducting  layer.  Reducing  the  doping  makes 
ohmic  contact  more  difficult  to  achieve,  while  thinning  the 
conducting  region  becomes  impractical  and  fringe  field 
effects  begin  to  contribute  as  the  substrate  p-n  junction 
becomes  sizeable  compared  to  the  surface  conducting  zone. 

All  of  these  measurements  require  a homogeneous 
interfacial  layer  to  produce  meaningful  results;42  sintered 
Al-Si  contacts  have  been  shown  to  exhibit  uneven  reaction 
near  the  edges  of  the  contact  with  the  formation  of 
excessive  eutectic  reaction  product,  and  contact  resistance 
measurements  with  various  sizes  of  contact  bar  show 
nonlinear  properties  consistent  with  this  edge  effect.  If 
uniform  contacts  are  not  achieved,  an  accurate  estimate  of 
the  interfacial  contact  resistance  cannot  be  made. 

Summary 

In  semiconductor  device  applications,  silicides  offer 
many  candidates  with  low  and  metal-like  resistivities  and 
high  temperature  stability.  They  are  also  expected  to  have 
good  resistance  to  electromigration  due  to  the  strong 
bonding  compared  to  metals,  and  if  formed  on  polysilicon 
they  can  be  useful  in  providing  lower  resistance  within  the 
current  MOS  gate  technology. 
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The  material  properties  of  interest  are  many,  but  of 
greatest  importance  to  device  design  are  sheet  resistance, 
Schottky  barrier  height  and  contact  resistance.  For 
stability  in  processing  and  device  reliability,  films  must 
be  stable  with  respect  to  the  substrate  and  dielectric 
layers  under  thermal  cycling  and  have  tolerable  oxidization 
behavior.  Physical  properties  of  films  such  as  stress, 
smoothness  and  adherence,  and  etchability  for  patterning  are 
also  of  concern.  The  properties  of  transition  metal 
silicides  are  discussed  below. 

Transition  Metal  Silicides 

Two  classes  of  transition  metal  silicides  have  emerged 
as  the  most  interesting  for  semiconductor  applications:  the 
refractory  metal  silicides  (metals  from  groups  IVB  to  VIB  in 
the  periodic  chart)  and  the  silicides  of  the  noble  and  near- 
noble metals  (group  VIII) . 

The  near-noble  and  noble  metals  are  so  named  because 
they  exhibit  little  or  no  affinity  for  oxygen,  while  the 
refractory  metals  have  high  melting  temperatures  and  are 
thus  often  suited  for  refractory  (high  temperature) 
applications.  The  two  classes  of  silicides  are  clearly 
differentiated  in  terms  of  the  temperature,  mechanism  and 
dominant  diffusing  species  in  the  thin  film  reaction  on 
silicon  substrates. 

They  also  display  characteristic  differences  in  terms 
of  the  properties  of  the  silicide  formed  and  its  interfacial 
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electronic  properties  with  respect  to  silicon.  While  the 
similarities  within  each  class  are  thus  emphasized,  it  is  to 
be  noted  that  there  are  exceptions,  as  well  as  properties 
that  vary  significantly,  and  not  always  in  a well-behaved 
fashion,  within  each  class  or  column  of  the  periodic  table. 

Formation  and  Reaction  Kinetics 

For  metal  films  on  silicon  or  polysilicon  substrates, 
the  metal-rich  phases  are  often  formed  first,  frequently 
consuming  all  of  the  metal  before  the  appearance  of  a phase 
richer  in  silicon.5  The  final  phase  is  generally  near 
M:Si::l:2  at  higher  temperatures  and  in  the  presence  of 
excess  silicon,  which  is  typical  of  the  processing 
environment  for  silicon  devices.  This  superior  stability  of 
the  dominant  phase  once  formed  on  silicon  is  advantageous 
from  the  point  of  view  of  device  processing. 

Studies  of  the  time  dependence  of  silicide  growth  in 
most  cases  indicate  a t1/2  dependence,  indicating  a 
diffusion-controlled  reaction.14  In  other  cases,  sometimes 
in  the  same  system,  a linear  time  dependence  was  observed, 
indicating  an  interface  rate-controlled  reaction,  raising 
the  questions  of  surface  cleanliness  and  clearly  leaving  the 
question  open  for  further  study. 

In  any  case,  the  refractory  metals  react  at  elevated 
temperatures  (T>600°C)  , and  silicon  is  the  dominant 
diffusing  species,  while  for  the  near-noble  metals  the 
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reaction  temperature  is  lower  (T<600°C)  , but  the  metal  atoms 
are  found  to  dominate  the  diffusion  and  perhaps  promote  the 
breaking  of  the  Si-Si  bonds  at  the  low  temperatures  due  to 
their  interstitial  migration.45 

The  tendency  for  the  reactions  to  proceed  to  the 
formation  of  silicon-rich  phases  can  be  understood 
thermodynamically,  in  that  the  heat  of  formation  per  metal 
atom  increases  with  increasing  silicon  content,  with  a 
disilicide  forming  a two-phase  field  with  silicon  at  the 
silicon-rich  end  of  the  binary  phase  diagram.14  However, 
efforts  to  formulate  rules  accurately  predicting  first-phase 
formation  and  second-phase  nucleation  from  phase  diagram 
parameters  have  not  been  totally  successful.  Factors 
affecting  first  phase  formation  include  cleanliness  of  the 
interface,  purity  of  materials,  diffusivity,  relative  free 
energy  of  formation  and  temperature  of  interaction.  It 
should  be  noted  that  the  presence  of  an  oxide  on  the 
interface  can  have  serious  effects  on  the  reaction,  even  for 
metals  with  low  affinity  for  oxygen.14 

For  the  refractory  metal  silicides,  the  crystal 
structures  of  the  final  silicon-rich  phases  follow  a trend 
across  each  period  in  the  chart,5  from  orthorhombic  (group 
IV)  to  hexagonal  (group  V)  to  tetragonal  (group  VI) . The 
near  noble  metals  are  more  disparate,  from  cubic  for  the 
third  period  (Co,  Ni)  to  hexagonal  (Pd)  and  orthorhombic 
(Pt)  . 
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Codeposited  silicides  are  generally  amorphous  as 
deposited.  On  silicon  or  polysilicon  they  nucleate  in  the 
disilicide  phase,  with  rejection  of  excess  silicon  or 
reaction  with  the  substrate  to  form  single-phase  disilicide. 
On  oxide,  they  tend  to  form  a phase  mixture  according  to  the 
composition  deposited.46 

Silicide  Properties 

Physical  properties 

Films  produced  by  the  thin  film  reaction  are 
potentially  subject  to  a large  tensile  stress,  due  to  the 
large  net  volume  shrinkage  resulting  from  silicide 
formation.14  The  relatively  open  silicon  lattice  is  capable 
of  accommodating  a large  portion  of  the  volume  of  metal 
atoms  introduced  in  the  formation  of  the  silicide.  Stress 
can  threaten  the  mechanical  stability  and  even  in  cases  the 
adhesion  of  films. 

However,  while  the  films  do  show  tensile  stress  at  room 
temperature  after  high  temperature  annealing,  it  is  much 
less  than  would  be  expected  from  the  volumetric  change  of 
the  reaction.  It  is  comparable  to  the  stress  level  of 
annealed  codeposited  stoichiometric  films,14  lxlO10  to  2xl010 
dyn/cm2. 

The  stress  in  codeposited  films  often  shows  a marked 
dependence  on  the  Si/metal  ratio.47  Codeposited  films  of 
correct  or  silicon-rich  stoichiometries  tend  to  have  low 
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stress  as  deposited,  but  annealing  produces  tensile  stresses 
in  these  films  as  well,  apparently  due  to  differences  in 
thermal-expansion  coefficients  with  respect  to  silicon 
substrates.48 

The  details  of  stress  development  vary  among  these 
metals.  Ti  and  Ta  films  are  compressive  as  deposited,  but 
with  annealing  develop  tensile  stress  which  levels  off  at 
high  temperatures.49  Co  films  are  in  tension  as  deposited, 
and  with  reaction  this  tension  increases,  then  is  reduced 
after  peaking  at  900°C5.  For  Pt  and  Pd  films  the  stress  is 
low  and  tensile  as  deposited,  increasing  with  annealing 
temperature,  and  leveling  off  or  decreasing  above  700°C.5Q 
Ni  films  have  a very  low  tensile  stress  that  increases 
slightly  with  temperature.  This  low  stress  may  be  related 
to  the  growth  of  epitaxial  NiSi2  on  silicon,  which  is 
observed  in  these  films.51,52 

For  TaSi2,  TiSi2,  MoSi2  and  WSi2,  the  silicides  show  no 
hysteresis  in  stress  as  a function  of  temperature  in  situ  in 
both  the  heating  and  cooling  cycles,  which  indicates 
stability  once  formed.48 

The  silicides  show  good  adhesion  for  the  stable 
silicide.  However,  at  temperatures  of  500°  and  600°C,  Ni 
films  and  Co  films  (thicker  than  lOOnm)  form  a monosilicide 
which  lifts  off  the  substrate.5 
Etchability 

Etching  of  silicides  is  problematical.  Wet  chemical 
etching  is  generally  ineffective  except  in  HF,  and  the  HF 
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will  attack  silicon  and  Si02  structures  in  the  device  as 
well.  In  the  lack  of  selective  chemical  etchants,  the 
choices  are  narrowed  to  1)  chemically  etching  a pattern  into 
a deposited  metal  before  annealing  to  form  the  silicide,  2) 
lift-off  of  a patterned  resist,  or  3)  reactive  ion  etching 
(RIE) . Developments  in  the  latter  have  the  most  long-range 
promise. 14 
Oxidation  behavior 

Device  processing  often  involves  oxidation  steps  to 
grow  an  insulating  layer  of  Si02  on  the  exposed  silicon. 

The  silicides  should  survive  these  oxidation  treatments,  and 
hopefully  would  be  capable  of  producing  a dielectric  layer 
on  the  silicide  as  well.  Lastly,  metallizations  must  be 
stable  on  Si02  layers  during  heat  treatments. 

Oxides  of  the  metals  show  melting  points  that  decrease 
from  group  IVB  to  VIB,  and  increase  with  each  period.5  Thus 
of  group  VIB  only  W03  has  a melting  temperature  above 
1000°C,  and  only  those  in  group  VB  heavier  than  V meet  this 
test. 

Volumetric  change  in  the  metals  upon  oxidation  also  do 
not  favor  group  VIB,  with  the  metals  exhibiting  a much 
larger  percent  volume  gain  compared  with  Si02,  which  would 
produce  large  accommodating  stresses  upon  oxidation.  Groups 
IVB  falls  below  Si02,  and  group  VB  slightly  above,  with 
respect  to  expansion  in  oxidation.5 

A comparison  of  the  heats  of  formation  of  metal  oxides 
versus  that  for  Si02  shows  that  for  group  VIB  metals,  the 
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oxides  are  less  stable  than  Si02  and  therefore  oxides  of  Cr, 
Mo  and  W can  be  reduced  by  Si.  Under  oxidation  in  the 
presence  of  excess  silicon,  silicides  of  these  metals  will 
be  stable;  when  there  is  no  excess  silicon  available,  they 
will  tend  to  lose  silicon  to  the  oxide  and  become  metal- 
rich.5 

Group  IVB  metals  can  reduce  Si02,  and  to  a lesser 
degree  so  can  the  group  VB.  Silicides  of  these  elements  may 
form  both  oxides.  For  other  transition  metals,  as  for  the 
group  IVB,  oxidation  of  the  metal  in  preference  to  the 
silicon  can  lead  to  decomposition  of  the  silicide,  with  a 
buildup  of  silicon  residue.5 

Experimental  results  in  general  concur  with  this 
thermodynamic  picture.  An  Si02  cap  grows  on  silicides  that 
have  an  intervening  layer  of  polysilicon  between  the 
silicide  and  an  underlying  oxide,  and  the  polysilicon  is 
consumed  without  affecting  the  silicide,  except  for  MoSi02 
and  WSi02  films,  which  as  expected  become  metal-rich  under 
the  Si02  cap.14  For  silicides  directly  on  the  oxide,  TiSi2 
and  TaSi2  are  reasonably  stable;  again  MoSi02  and  WSi02  are 
unstable,  producing  metal-rich  films,  metal  in  the  films, 
and  evolving  volatile  metal  oxides.  For  gate  electrodes 
only  silicides  such  as  those  of  Ta  and  Ti  seem  recommended. 
Sheet  resistance 

A great  deal  of  the  interest  in  silicides  is  due  to 
their  low  resistivity  with  respect  to  conventional 
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polysilicon  gate  metallizations.  Silicides  of  the 
transition  elements  vary  from  near  15nQcm.  to  over  lOOOnQcm. 
in  resistivity.14  For  thin  films,  resistivities  are 
generally  higher  than  bulk  values,  but  a number  of  reacted 
silicides  have  been  reported  in  the  promising  range  of  40- 
60nOcm.  Codeposited  films  are  generally  more  resistive  than 
reacted  films  of  the  same  element. 

For  the  silicides  formed  from  metals  of  group  IVB  to 
VIIB  of  the  periodic  table,  it  has  been  observed  that  in  a 
given  period  an  increase  in  resistivity  is  observed  with 
increasing  number  of  d-electrons.  This  change  is  associated 
with  a change  in  the  crystal  structure,  from  orthogonal  to 
hexagonal  to  tetragonal,  and  both  changes  are  attributable 
to  the  changing  electronic  structure  of  the  metal  atom.  It 
has  been  speculated  that  the  very  low  resistivity  of  the 
group  IVB  metal  silicides  is  related  to  a hybridization  of 
the  two  d-electrons  with  s-p  electrons  from  the  silicon. 53,54 
With  increasing  metal  atom  size,  the  periodic  change  is  less 
pronounced,  and  for  group  VIII  elements  the  resistivity  is 
low  and  relatively  constant. 

By  comparison,  the  resistivity  of  the  transition  metals 
themselves  tends  to  decrease  across  each  row,  with  an 
anomalous  peak  at  group  VIIB  attributed  to  the  very  stable 
d5  configuration,  and  group  VIII  metals,  like  the  silicides, 
have  low  resistances  that  are  similar  in  value. 


34 


Rare  Earth  Silicides 

The  thin  film  reaction  of  a silicon  substrate  with  a 
deposited  rare  earth  metal  produces  a stable  and  homogeneous 
disilicide  film  with  a very  low  Schottky  barrier  height  to 
n-silicon6'7  (0.3-0. 4 eV)  [Table  2-1].  These  low  barriers 
suggest  that  low-resistance  ohmic  contacts  may  be  achieved 
with  electron  tunneling  as  a dominant  conduction  mechanism. 
Given  their  high  melting  temperatures55,56  [Table  2-2]  the 
rare  earth  silicides  are  clearly  excellent  candidates  for 
low  resistivity,  thermally  stable  contacts  for  n-Si. 

Thin  Film  Reaction  on  Silicon 

The  rate  of  the  rare  earth-silicon  reaction  is 
difficult  to  determine  due  to  the  extreme  rapidity  with 
which  the  metal  is  consumed.45  Using  ion-beam-induced 
mixing57  or  rapid  electron-beam  pulses  to  dissolve 
interfacial  impurities58  permitted  measurement  of  the 
activation  energy  for  erbium-silicon  (Ea=1.75  eV)  and  Tb-Si 
(Eg=l . 6 eV)  reactions.  Thus,  even  at  the  very  low 
temperatures  at  which  the  reaction  is  initiated  through  the 
interfacial  barrier,  rates  of  formation  of  the  silicide  are 
rapid.  The  result  is  the  rough  and  pitted  surface  that  is 
typical  of  the  rare  earth/silicon  thin  film  reaction. 9,59,60 

In  addition,  it  has  been  pointed  out  that  rare  earth 
disilicides  can  form  a hexagonal  phase  (A1B2  type)  with  a 
[0001]  face  that  is  closely  related  to  the  (111)  face  of 
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silicon.10  Epitaxial  thin  films  have  been  demonstrated  for 
the  higher  lanthanides  (Z>63)  using  rapid  anneals  by 
electron  beam  or  laser.11'12  Epitaxy  of  the  orthorhombic 
"GdSi2"  allomorph  on  Si  (100)  has  also  been  reported.13  Such 
epitaxial  layers  may  offer  potential  advantages  in  interface 
stability  and  properties,  for  example  enhanced  resistance  to 
diffusion  (i.e.  diffusion  barriers)  if  penetration  of  the 
silicide  by  conductive  overlayers  is  reduced  by  the 
elimination  of  grain  boundaries. 

Rare  Earth  Contacts  to  GaAs 

Rare  earth  silicide  contacts  to  GaAs  have  not  been 
reported.  However,  photoemission  studies  of  several  rare 
earth  metals  on  cleaved  GaAs ( 110) 61"65  show  room  temperature 
reaction  resulting  in  an  abrupt  interface  and  Fermi-level 
pinning  in  the  bandgap  [Table  2-3]. 

The  nature  of  Schottky  barrier  formation  continues  to 
be  controversial,66  with  schools  supporting  the  roles  of 
defect  states,67,68  chemical  reactivity  and  effective  work 
function  of  interfacial  compounds,69  and  strength  of  the 
chemical  bond  and  its  dipole70  as  factors  in  the  pinning 
behavior.  The  high  reactivity  and  low  electronegativity  of 
the  rare  earth  metals  makes  them  good  test  cases  for  these 
divergent  theories,  but  the  barrier  heights  measured  do  not 
show  agreement  beyond  the  expected  band-gap  pinning. 

In  so  far  as  the  barrier  height  of  the  silicide  will 
follow  that  of  the  metal,  and  given  the  inevitable  pinning 
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in  the  band  gap  for  GaAs,  successful  Ohmic  contacts  will 
depend  upon  tunneling  or  thermionic  field  emission  as 
dominant  conduction  mechanisms.  Calculations  have  shown 
that  these  mechanisms  can  produce  low  contact  resistances  on 
highly  doped  (nD^1018  cm'3)  substrates,30  with  exponential 
decrease  in  Rc  with  decreasing  barrier  height.  Thus  the 
discrepancy  in  barrier  heights  reported  for  rare  earth 
metals  [see  Table  2-3]  would  translate  to  a factor  of  5 to 
10  in  contact  resistance  for  tunneling  Ohmic  contacts.  This 
underscores  the  significance  of  these  barrier  height 
measurements  for  Ohmic  contacts. 

Resistivities 

Bulk  and  thin  film  measurements  show  moderate 
resistivities  for  the  rare  earth  silicides9,56,71  [Table  2-4, 
2-5],  which  indicates  that  they  will  be  better  suited  to 
multilevel  metallization  schemes  with  a more  conductive 
metallization  layer.  Interestingly,  semiconducting 
properties  and  a band  gap  of  0.19  eV  have  been  reported  for 
LaSi2. 72,73  If  this  were  true  of  other  silicides  it  would 
have  implications  for  performance  in  devices. 

Problems  with  Reacted  Thin  Films 

The  main  problems  encountered  in  the  previous  studies 
of  reacted  thin  films  on  silicon  are  related  to  oxygen 
contamination,  either  incorporated  oxygen  in  the  deposited 
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rare  earth  metal  or  surface  oxides  on  the  silicon  substrate, 
which  are  dissolved  by  the  strongly-reducing  rare  earth 
metal.  High  concentrations  of  oxygen  in  the  metal  inhibit 
the  onset  of  the  silicide  reaction,  and  have  been  also  been 
identified  as  the  cause  of  the  extreme  roughness  of  the 
resulting  films,  due  to  rapid  reaction  beginning  at  points 
of  weakness  in  an  oxide  barrier  layer.74 

This  pitting  has  been  identified  as  a factor  in  poor 
diode  behavior  and  is  undesirable  for  the  fabrication  of 
shallow  contacts.75  The  pitting  has  been  reduced  or 
eliminated  by  ion-implantation  or  mixing,57  rapid  annealing 
by  electron  beam76  or  halogen  lamp,12  codeposition  of  metal- 
silicon  mixtures,77  UHV  deposition  on  atomically  clean 
surfaces,78  or  by  prereaction  of  very  thin  metal  films 
before  depositing  the  desired  thickness  of  rare  earth 
metal.79 

Crystallographic  Structures  of  Rare  Earth  Silicides 

Rare  earth  silicides  near  the  disilicide  composition 
may  exhibit  several  crystal  structures,  depending  upon 
temperature,  composition,  size  of  the  metal  atom,  and 
valency  of  the  compound.80  The  complexity  of  rare  earth 
crystallography  makes  phase  identification  difficult  for 
these  films.59 

Gadolinium  with  its  characteristic  orthorhombic  "GdSi2" 
structure  lies  at  the  crossing  point  of  tendencies  to  form 
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tetragonal  ThSi2-type  structures  (lower  lanthanides)  or  to 
form  hexagonal  structures  of  AlB2-type  or  a silicon- 
deficient  "defect"  A1B2  (higher  lanthanides)  .81 

The  hexagonal  lattice  can  be  arrived  at  by  a small 
distance  shift  of  a few  atoms  along  the  [010]  direction  of 
the  tetragonal  lattice82  [Figure  2-1],  which  can  explain  the 
reversible  temperature-dependent  transformation  observed  in 
these  compounds.83,84*85  The  hexagonal  form  is  favored  at 
lower  temperatures,  low  silicon  content,  and  low  population 
of  free  electrons  (low  valency) , as  well  as  low  atomic  size 
(higher  atomic  number) . The  orthorhombic  "GdSi2"  distortion 
of  the  tetragonal  lattice  is  an  intermediate  form,  also 
characterized  by  a deficiency  in  silicon  and  intermediate 
atomic  size  ratio  (rS)-/rM  = 0.58). 81  At  higher  temperatures 
the  tetragonal  ThSi2-type  structure  has  also  been  reported 
for  Gd  and  other  rare  earths82,86,87  [Table  2-6]. 

Both  the  orthorhombic  and  the  hexagonal  forms  have  been 
observed  in  reacted  films  of  rare  earth  silicides,59 
although  Gd  favors  the  orthorhombic  form  both  in  bulk83  and 
thin  film  samples.9  An  exception  is  found  for  liquid-phase 
reaction  of  Gd  films  on  Si  (111)  by  electron-beam  melting. 11,12 
Hexagonal  GdSi,  7 has  been  also  been  formed  in  bulk  samples 
using  a low-temperature  amalgam  method  (T<600°C)81  or  by 
higher  temperature  reaction  of  silicon-deficient  mixtures.87 

The  hexagonal  form  of  disilicide  exhibits  two  distinct 
AlB2-type  structures  for  the  actinide  and  lanthanide  metals, 
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one  stoichiometric  and  the  other  deficient  in  silicon.82 
For  the  silicon  atoms  to  fit  perfectly  into  interstitial 
sites  in  a hexagonal  A1B2  structure  requires  a low  size 
ratio  of  the  metal  atom  to  the  silicon  atom  (rsi/rM  ^ 0.527). 
Since  these  compounds  that  do  not  satisfy  the  size  ratio 
requirement  for  ideal  interstitial  configuration,  they  form 
a defect  A1B2  structure  in  which  the  unit  hexagonal  cell 
must  expand  in  either  the  [0001]  or  the  [1000]  direction. 
Stoichiometric  compounds  follow  the  latter  arrangement,  and 
as  a result  the  axial  ratio  of  the  cell  c/a  is  less  than 
unity  (i.e.,  for  ideal  interstitial  A1B2,  c/a  = 1).  For 
nonstoichiometric  or  defect  A1B2  structures,  the  expansion 
is  in  the  axial  direction,  so  that  c/a  > 1. 

The  rare  earth  silicides  form  the  AlB2-type  structure 
only  when  the  atomic  size  ratio  rsj/rM  is  greater  than  0.579, 
well  above  the  required  ratio  for  interstitial  arrangement 
(rS)-/rM  ^ 0.527).  These  compounds  exhibit  an  axial  expansion 
of  the  unit  cell,  and  a tendency  to  be  nonstoichiometric. 
Thus  Gd3Si5,  in  the  hexagonal  form,  has  axial  ratio  c/a  = 
1.076. 87  Since  disilicides  of  Tm  and  Yb  have  been  observed 
to  maintain  this  structure  even  at  2:1  stoichiometry,84  it 
appears  that  the  alternative  radially  expanded  A1B2  lattice 
does  not  occur  for  the  rare  earth  silicides. 

For  the  orthorhombic  "GdSi2"  form,  considerable 
deviation  from  stoichiometry  has  been  observed  in  bulk 
samples.  Perri  et  al.86  found  that  both  density 
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measurements  and  X-ray  intensity  analysis  indicated  that  the 
stoichiometry  was  GdSi., 4,  and  lower  density  has  also  been 
found56  for  the  orthorhombic  DySi2  [Table  2-7].  More  recent 
RBS  studies  of  orthorhombic-phase  GdSix  thin  films  indicate 
stoichiometry  from  1.4:1  to  2:1,  a considerable  range  of 
homogeneity,  with  stable  properties  beginning  at  GdSi1  7.77 
There  remains  potential  confusion  concerning  this  point 
because  many  studies  seem  to  rely  on  measurements  of  the 
hexagonal  A1B2  phase  which  is  also  commonly  assigned  a ratio 
of  1.7:1. 88  In  this  study,  reacted  products  of 
stoichiometry  1.7:1  were  chosen  as  the  baseline  for 
composition. 
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Table  2-1  Schottky  barrier  heights  [in  eV]  of  several  rare 
earth  disilicides  to  silicon. 


YSi2 

GdSi2 

oysi2 

HoSi2 

ErSi2 

ref. 

n-Si 

0.39 

0.37 

0.37 

0.37 

0.39 

6 

p-Si 

0.75 

0.71 

0.73 

— 

0.70 

0.77 

6 

7 

Table  2-2  Melting  temperatures,  GdSi2. 


Compound 

T^c° 

ref. 

GdSi2 

1500 

55 

GdSi2 

1540 

56 

Table  2-3  Schottky  barrier  heights  [in  eV]  of  several  rare 
earth  metals  to  GaAs. 


Metal 

ref. 

Sm 

0.63 

64 

Tb 

0.85 

61 

Dy 

0.85 

61 

Er 

0.85 

61 

Yb 

0.61- 

0.72 

65 
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Table  2-4  Bulk  resistivities  for  rare  earth  silicides. 


K 

Silicide 

[(iQcm.  ] 

ref. 

ScSii.6 

40.6 

71 

YSi2 

69.0 

71 

LaSi2 

350.0 

71 

273 

56 

CeSi2 

183.0 

71 

NdSi2 

349 

56 

SmSi2 

210.0 

71 

GdSi2 

263 

56 

Table  2-5  Thin  film  resistivities  of  rare  earth  silicides  and 
rare  earth  metals,  and  bulk  resistivities  of  the  metals,  based 
on  Table  I.,  ref.  9. 


Metal 

Y 

La 

Gd 

Dy 

Ho 

Er 

Thin  film 
R0 , jiDcm . 

95 

104 

133 

116 

99 

81 

Bulk 
R0 , ^cm . 

60.2 

79.4 

135.9 

92.6 

80.6 

85.5 

Silicide 

YSi2 

LaSi2 

GdSi2 

DySi2 

HoSi2 

ErSi2 

R0 , jiGcm . 

135 

116 

207 

204 

185 

115 
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Table  2-6  Unit  cell  dimensions  for  the  three  reported 
structures  of  GdSix  (1.7^x^2). 


Crystal  type 

a 

b 

C 

ref . 

Tetragonal 

4.10 

13 . 61 

83 

Orthogonal 

4.09 

4.080 

4.01 

3.996 

13.44 

13.41 

86 

80 

Hexagonal 

3.877 

3.864 

4.172 
4 . 155 

87 

80 

Table  2-7  Anomalous  densities  for  orthorhombic  rare  earth 
silicides. 


Silicide 

Density 
( ideal) 
[g/cm3] 

ref. 

Density 
(meas . ) 
[g/cm3] 

ref. 

DySi2 

6.8 

55 

5.2 

56 

GdSi2 

6.4 

55 

5.94 

86 
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Figure  2-1  Hexagonal  and  tetragonal  lattices  for  ThSi2 

(a)  S-ThSi2  (hexagonal  AlB2-type)  packing  diagram. 

(b)  Tetragonal  structure  cell  of  cr-ThSi2  projected  along 
[010]  direction. 

(c)  Hexagonal  structure  cell  of  £-ThSi2  projected  along 
[0001]  direction. 


CHAPTER  THREE 

DEPOSITION  AND  ANALYTICAL  TECHNIQUES 
Introduction 

Ion-beam  sputter  deposition  from  a composite  target  of 
silicon  and  gadolinium  was  used  to  produce  GdSix  films  on 
Si,  GaAs  and  amorphous  quartz  slides  to  study  the  properties 
of  codeposited  mixtures.  For  comparison,  sample  thin  films 
of  gadolinium  disilicide  were  produced  by  thin  film  reaction 
of  electron-beam  deposited  rare  earth  metal  with  substrate 
silicon  or  deposited  amorphous  silicon,  but  the  times, 
temperature,  and  completion  of  the  reaction  were  so 
unpredictable  because  of  oxygen  contamination  that  this 
approach  was  found  impractical  for  producing  silicides. 

Compositions  of  the  codeposited  films  were  varied  from 
N„./N„,  of  14:1  to  1:1.  These  films  were  then  studied  in 
comparison  to  the  thin  film  reaction  product  and  to  a bulk- 
prepared  disilicide  from  Cerac,  Inc. 

The  composition  was  evaluated  using  Rutherford 
backscattering  spectroscopy  (RBS)  and  Auger  electron 
spectroscopy  (AES) .89  The  AES  instrument  used,  having 
enhanced  imaging  and  areal  resolution  capabilities,  is  also 
known  as  a scanning  Auger  electron  microprobe  (SAM) . 
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Post-analysis  using  SAM  showed  that  divergent  RBS  data 
resulted  from  topography  in  the  reacted-film  and 
heterogeneous  phase  formation  in  the  bulk-prepared  samples, 
respectively.  AES  spectra  from  the  homogeneous  regions  of 
the  two  silicides  was  used  to  evaluate  homogeneous  Gd/Si 
films  deposited  from  the  composite  sputter  target. 

Evaluation  of  surface  composition  gradients  and  the 
development  of  quantitative  Auger  electron  spectroscopy  for 
the  Gd/Si  system  will  be  presented  separately,  in  Chapters 
Four  and  Five. 

The  structure  of  the  films  were  studied  using  x-ray 
diffraction  ( XRD ) 89  and  also  electron  diffraction  in  the 
transmission  electron  microscope  (TEM) . Other  measurements 
obtained  included  the  resistivity  (by  the  four-point  probe 
technique) , optical  properties  (using  optical  absorption  and 
reflectometry) , and  contact  properties  (using  current- 
voltage  measurements  of  0 . 7mm-diameter  dot  contact  diodes, 
deposited  through  a mechanical  mask) . These  results  are 
presented  in  Chapter  Six. 

Thin  Film  Deposition 

All  substrates  were  degreased  using  a standard  cleaning 
procedure  (TCE  - Acetone  - Methanol) . Silicon  and  GaAs 
substrates  were  then  chemically  etched  in  a solution  of  1:10 
HF:H20  to  remove  surface  oxides.  Samples  for  ion-beam 
deposition  were  further  cleaned  under  vacuum  in  the 
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deposition  chamber  by  exposing  the  sample  surface  to  a flux 
of  =300eV  Ar+  ions  before  deposition. 

Ion-beam  sputter  codeposition 

The  problems  inherent  in  the  thin  film  reaction 
approach  provided  incentive  to  develop  a codeposition 
process  in  which  the  silicon  and  gadolinium  are  deposited  as 
a stoichiometric  mixture,  thus  avoiding  both  excessive 
oxidation  of  the  metal  before  heat  treatment  and 
irregularities  in  the  films  due  to  uneven  redistribution  of 
species  in  the  course  of  the  reaction.  Besides,  the 
cosputtering  technique  was  desirable  in  terms  of 
demonstrating  reasonable  film  qualities  that  could  be 
transferred  to  applications  without  elaborate  techniques 
that  are  time-consuming  and  expensive.  Finally,  the 
cosputtering  approach  was  made  more  attractive  by  the 
development  of  a target  configuration  that  allowed  for 
tuning  the  Gd/Si  ratio  by  altering  the  orientation  of  the 
composite  target  of  Si  and  Gd  [Figure  3-1].  Compositions  of 
the  films  were  varied  from  Nsj/NGd  of  14:1  to  0.87:1,  and  by 
comparative  analysis  the  deposition  of  films  near  the 
desired  GdSi,  7 composition  was  maintained.  The  codeposited 
films  were  free  of  the  roughness  and  pitting  typical  of  the 
reacted  films. 

The  deposition  chamber  was  a conventional  VEECO 
Microetch  system  consisting  of  a diffusion-pumped  glass  bell 
jar,  with  a base  pressure  of  =lxl0'6  Torr,  backfilled  with 
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5xl0"6  Torr  of  Ar  during  sputtering.  The  ion  source 
produces  a uniform  parallel  beam  of  ions  with  an 
accelerating  voltage  adjustable  in  the  0-3keV  range, 
allowing  substrate  cleaning  with  lower-energy  ions  while 
depositions  were  done  at  2keV.  A quartz  crystal  thickness 
monitor  mounted  in  the  chamber  was  calibrated  and  cross- 
checked by  profilometer  measurements  on  the  deposited  films. 
Mechanical  shadow  masks  fabricated  from  tantalum  metal  foil 
were  used  to  define  circular  diode  test  patterns  as  required 
for  contact  measurements. 

Silicides  formed  by  reaction  of  deposited  Gd  metal 

In  order  to  assess  the  properties  of  silicides  produced 
by  the  thin  film  reaction  approach,  thin  films  of  Gd  metal 
were  prepared  by  electron  beam  evaporation  and  magnetron 
sputter  deposition  for  reaction  studies,  both  on  silicon  and 
GaAs  substrates. 

Films  of  Gd  metal  were  deposited  by  electron  beam 
evaporation  in  an  ion-pumped  chamber,  with  turbo-pumping  to 
assist  in  handling  the  gas  load  during  evaporation.  The  base 
pressure  in  the  system  was  =lxl0'6  Torr  and  the  pressure 
during  evaporation  was  ®5xl0'5  Torr. 

Samples  were  coated  to  various  thicknesses  with  Gd  or 
Gd/Si  layers.  All  samples  were  capped  with  =20nm  of  silicon 
to  protect  the  surface  from  oxidation  during  sample  transfer 
to  analysis. 
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Additional  samples  of  Gd  on  Si  were  prepared  by  DC 
magnetron  sputter  deposition  in  a diffusion-pumped  chamber, 
base  pressure  =lxl0'6  Torr,  backfilled  with  Ar  gas  during 
deposition.  The  substrates  were  chemically  cleaned,  and  then 
plasma  cleaned  in  the  chamber.  Again,  they  were  capped  with 
silicon  (=20nm)  before  removing  them  from  the  chamber. 

The  thin  film  reaction  has  several  disadvantages  as  a 
technique  for  producing  silicide  films.  The  deposited  rare 
earth  metal  films  are  extremely  susceptible  to  oxidation7,59 
and  the  reaction  has  a tendency  to  produce  an  irregular  and 
pitted  interface  with  the  silicon.  Good  results  from  the 
reaction  require  careful  attention  to  the  cleanliness  of  the 
substrate  surface  and  annealing  under  vacuum  or  a strongly 
reducing  atmosphere. 

For  measurement  of  some  intrinsic  properties  and  for 
applications  to  compound  semiconductors,  it  is  desirable  to 
produce  the  silicide  on  substrates  other  than  silicon.  This 
would  require  alternate  or  multiple  layers  of  Si  and  rare 
earth  metal  of  appropriate  thickness  ratio  with  subsequent 
reaction  to  form  the  silicide.  The  reaction  of  deposited 
gadolinium  with  overlayers  of  amorphous  silicon  does  offer 
advantages  in  that  reaction  has  been  shown  to  occur  at  lower 
temperatures  and  to  produce  smoother  films  than  the  reaction 

r 

with  crystalline  substrates. 

Some  attempts  were  made  to  create  reacted  films  on  non- 
silicon substrates  but  were  abandoned  since  the  difficulty 
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of  producing  uniform  films  of  correct  stoichiometry  using 
this  approach  was  difficult,  and  optical  analysis  was  a 
priority  that  required  this  uniformity.  This  approach 
requires  the  deposition  of  a bilayer  (or  multilayer)  of  Gd 
and  Si  in  appropriate  thickness  ratio,  or  alternatively 
deposition  of  a surplus  of  Si  as  a cap  for  the  annealing 
process  (subsequently  to  be  removed) . 

However,  the  technique  is  unattractive  in  terms  of 
controlling  processing  parameters.  That  is,  while  such 
multilayer  structures  have  the  potential  for  reducing  or 
eliminating  the  pitting  and  unevenness  of  the  films,  and 
could  produce  silicide  films  on  substrates  other  than 
silicon,  considerable  complications  are  introduced  in 
controlling  the  deposition  and  reaction. 

Analytical  Instruments 

RBS . 89  The  samples  were  analyzed  using  2.0  MeV  alpha 
particles  with  a beam  current  of  120  nAmp.  Spectra  were 
smoothed  and  stoichiometries  determined  by  fitting  with 
self-consistent  stopping-power  calculations.  The  chamber 
was  equipped  with  a resistive  heating  stage  and  reaction  of 
thin  films  of  Gd  on  Si  substrates  was  monitored  in  situ. 

SAM.89  The  AES  data  were  acquired  on  a PHI  660 
Scanning  Auger  Microprobe  in  the  N(E)  mode.  Base  pressure  in 
the  analytical  chamber  was  =4xlO'10  Torr,  with  primary  beam 
energy  of  lOkeV  and  a beam  current  =400  nAmp.  The  surfaces 
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were  sputter-cleaned  using  a differentially  pumped  ion  gun. 

- A 

Chamber  pressures  during  sputtering  were  =8x10  Torr. 

The  areal  resolution  and  imaging  capabilities  of  the  SAM 
made  it  very  useful  in  comparing  the  silicides  prepared  by 
the  various  techniques. 

XRD . 89  Measurements  were  carried  out  with  a Phillips 
powder  diffractometer  using  Cu  Kff  radiation  and  the 
crystallization  properties  of  the  films  were  studied  using  a 
hot-stage  sample  holder. 

Composition  of  Films 

The  reacted  thin  film  specimens  along  with  bulk- 
prepared  GdSi2  provided  a reference  for  comparison  with  the 
codeposited  Gd-Si  films.  Relative  heights  of  the  leading 
edges  of  the  Si  and  Gd  elemental  transitions  in  the  RBS 
spectra  was  used  as  an  alternative  method  for  determining 
the  compositions  [Figure  3-2].  Peak  height  ratios  (Gd/Si) 
and  integral  peak  areas  were  not  entirely  consistent  among 
the  three  types  of  samples. 

Subsequent  examination  of  microstructure  in  the  bulk 
disilicide  by  SAM  showed  platelets  of  a silicon-rich 
precipitate  which  has  been  identified  as  silicon  carbide. 

The  presence  of  this  second  phase  is  believed  to  result  from 
reacting  a silicon-rich  2:1  Si:Gd  mixture  in  a graphite 
crucible,  while  the  stable  silicide  is  GdSi,  7.80  When  the 
matrix  phase  of  the  bulk  disilicide  was  analyzed  (as 


52 


compared  to  the  excess  silicon  phase  in  the  granules)  by 
SAM,  it  compared  favorably  with  the  fully  reacted  thin  film 
silicides  and  with  films  of  1.7:1  (RBS)  stoichiometry 
deposited  by  sputtering. 

It  was  also  noted  that  samples  produced  by  thin  film 
reaction  exhibited  a very  rough  surface  and  pitted  interface 
with  the  substrate.  Neither  of  these  fit  the  RBS  model  of  a 
planar  layerwise-homogeneous  sample.  The  codeposited  films, 
which  physically  resemble  the  RBS  model,  were  evaluated  for 
uniformity,  and  compositions  were  calculated  from  the  RBS 
model . 

Thus  multiple  phases  and  topography  effects  were  found 
to  cause  disagreement  between  the  SAM  and  RBS  for  the  bulk 
silicide  and  the  reacted  films.  For  the  smooth  homogenous 
films  deposited  by  the  ion  beam  technique,  RBS  and  SAM 
stoichiometry  data  can  be  compared  directly  and  showed  good 
agreement  after  application  of  matrix-corrected  AES 
sensitivity  factors,  which  are  developed  in  Chapter  Seven. 
The  corrected  AES  measurements  facilitated  the  monitoring  of 
composition  in  subsequent  depositions. 

The  polymorphism  of  the  rare  earth  disilicides 
mentioned  above  has  some  significance  to  the  questions  of 
composition  and  structural  analysis  for  these  films,  and 
will  be  discussed  below  in  Chapter  Six. 
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(b) 


Figure  3-1  Codeposition  from  a composite  sputter  target. 

(a)  Profile  view  of  diode  ion-beam  sputter  apparatus. 

(b)  Face-on  view  of  the  target,  showing  the  image  of  the  ion 
beam  on  the  target,  with  composition  of  the  sputtered  flux 
controlled  by  rotation  of  the  target. 
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Figure  3-2  RBS  analysis  of  gadolinium  disilicides. 

(a)  Thin  film  reaction:  1)  RBS  spectrum  of  a Gd  metal  film 
as  deposited  on  a Si  substrate,  and  2)  after  reaction 
(T=400°C) . In  formation  of  the  silicide,  the  silicon  peak 
has  shifted  to  the  surface  with  a lower  peak  height;  the 
gadolinium  peak  is  reduced  in  height  and  broadened. 

(b)  Codeposited  Gd-Si  mixture:  stoichiometry  is  determined 
by  comparing  the  integral  counts  in  the  silicon  shoulder  and 
in  the  Gd  peak.  These  samples  provide  the  most  accurate 
results  because  they  more  closely  conform  to  the  RBS  model 
of  planar,  layerwise  homogeneous  samples. 

(c)  Bulk-prepared  disilicide:  comparison  of  peak  heights  of 
the  leading  edges  of  the  Gd  and  Si  transitions  can  be 
compared  with  the  thin  film  specimens. 


CHAPTER  FOUR 

ANGLE-RESOLVED  X-RAY  PHOTOELECTRON  SPECTROSCOPY  (ARXPS) 

Introduction 

Formation  of  contacts  by  vacuum  deposition  requires  an 
understanding  of  the  condition  of  surfaces.  Chemically 
cleaned  substrates  quickly  develop  an  oxidized  surface  layer 
in  the  ambient.  Sputter  cleaning  can  change  the  surface 
composition  of  a compound  semiconductor  substrate  and 
preferential  sputtering  affects  the  results  of  film  analysis 
as  well. 

Electron  spectroscopies  produce  compositional  data  by 
measuring  the  relative  intensities  of  photoelectrons  or 
Auger  electrons  emitted  by  a sample  under  an  excitation  beam 
of  high  energy  photons  or  electrons.  The  signal  electrons 
counted  are  those  that  escape  the  sample  without  significant 
loss  of  energy  from  their  characteristic  excitation  peak. 

The  depth  analyzed  is  thus  related  to  the  inelastic  mean 
free  path,  that  is  the  average  distance  traveled  without  an 
energy-loss  scattering  event. 

Angle-Resolved  X-ray  Photoelectron  Spectroscopy  (ARXPS) 
has  been  recognized  as  a potential  method  for 
non-destructive  determination  of  a composition  depth  profile 
of  a surface  region  that  is  dimensionally  comparable  to  the 
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mean  free  path  of  the  photoelectrons.  Typical  values  for  the 
mean  free  path  are  on  the  order  of  10  to  50  A.  The  ability 
to  resolve  composition  changes  within  so  narrow  a surface 
region  combined  with  the  non-destructive  nature  of  the 
analysis  makes  the  technique  particularly  attractive  for 
studies  of  thin  films,  interfaces,  and  surfaces. 

An  inverse  Laplace  transform  method  has  been  used  to 
extract  composition  depth  profiles  from  angle-resolved  x-ray 
photoelectron  spectroscopy  (ARXPS)  data  from  oxidized  Si  and 
Ar+  bombarded  GaAs.  An  iterative  procedure  determines  a 
least-squares  curve-fit  of  the  data  using  the  Laplace 
transforms  of  step  functions,  which  are  summed  to  yield  the 
composition  profiles. 

In  this  chapter,  after  introducing  the  model  and  the 
iterative  approach  developed,  which  gives  a self-consistent 
solution  to  the  flux  equation,  two  example  applications  are 
presented.  1)  A sample  of  native  oxide  is  present  on  a Si 
wafer,  for  which  the  composition  profile  derived  defines  the 
thickness  and  chemical  layering  of  the  various  silicon  oxide 
states,  and  2)  samples  of  GaAs  have  been  sputtered  with  1.5, 
3.0,  and  5.0  keV  Ar+,  where  the  composition  depth  profiles 
show  a subsurface  As  depletion  the  extent  of  which  increases 
with  increasing  ion  energy,  resulting  from  preferential 
sputtering  caused  by  surface  segregation  enhanced  by 
sputter-assisted  diffusion  in  the  near  surface  region. 
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The  Inverse  Laplace  Transform:  An  Iterative  Solution 

A simplified  expression  of  the  photoelectron  flux 
generation  and  attenuation  in  the  sample  has  been  shown  to 
take  the  form  of  a Laplace  transform  of  the  depth 
composition  function.90  Unfortunately,  the  inverse  Laplace 
transform  relation  does  not  define  a unique  solution,  but 
rather  an  extensive  and  varied  set  of  solution  functions 
whose  transforms  are  nearly  identical.91  Consequently  the 
numerical  solution  of  the  composition  function  has  required 
elaborate  computational  techniques  combined  with  constraints 
on  the  form  and  the  range  of  the  solution  function. 

A new  approach  to  the  problem  has  been  introduced  in 
which  the  ARXPS  data  are  curve-fitted  using  the  concisely 
defined  Laplace  transforms  of  simple  step  functions.92 
Iterative  error  reduction  by  least-squares  comparison  yields 
a solution  consisting  of  summed  step  functions  as  an 
approximation  of  the  composition  function.  This  procedure 
does  not  involve  restrictions  on  the  form  of  the  solution 
and  reduces  considerably  the  complexity  of  the  computations 
involved.  In  this  work,  the  data  reduction  technique  has 
been  applied  to  ARXPS  data  from  two  different  types  of 
single  crystal  semiconductor  substrates  with  altered  surface 
layers.  First,  the  thin  native  oxide  on  a (100)  Si  substrate 
was  analyzed  using  the  Si  2p  photoelectron  peaks  of  the  XPS 
spectrum.  In  addition  to  the  core  line  representing  the 
elemental  Si  substrate,  several  shifted  peaks  were  resolved 
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with  binding  energies  identified  with  distinct  oxidation 
states  of  the  silicon  (i.e.,  SiO,  Si203,  and  Si02)  . The 
location  and  relative  abundance  of  Si  in  each  of  these 
oxidation  states  was  recovered  from  the  ARXPS  data.  The 
second  example  is  (100)  p-type  GaAs  (5xl018  cm'3  Zn)  analyzed 
after  1.5,  3.0,  and  5.0  keV  Ar+  ion  etching.  Changes  in 
surface  and  subsurface  composition  produced  by  the  ion 
bombardment  are  discussed  in  the  context  of  preferential 
sputtering. 

The  ARXPS  Flux  Model 

XPS  is  a surface  sensitive  analytical  technique  because 
of  the  rapid  energy  attenuation  of  low  energy  electrons  in  a 
solid.  This  surface  sensitivity  can  be  increased  by  changing 
the  take-off  direction  for  electrons  entering  the  analyzer 
toward  grazing  angles  relative  to  the  sample  surface  [Figure 
4-1].  This  decreases  the  effective  escape  depth  of  the 
photoelectrons,  quenching  the  signal  arising  from  deeper 
layers  of  the  sample,  thus  yielding  collected  data  that  are 
enriched  in  surface  information. 

The  intensity  of  photoelectron  flux  from  a given 
element  with  a concentration  Cj(x)  is  a function  of  the 
angle  d,  given  as93 

Ii  (6)  = ^(0)  Oi  f Cx(x)  exp  ( --  - * ) dx  (4-1) 

Jo  A jSUlD 
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For  a given  photoelectron  energy,  Ki  (9)  represents  a 
combination  of  the  effective  analyzer  response  function  and 
the  excitation  flux  as  a function  of  9.  The  product  of  this 
excitation  flux,  the  depth  composition  function  C^x),  and 
the  photoionization  cross-section  0-,  is  a generation 
factor.  The  exponential  term  represents  the  attenuation  of 
the  flux  contribution  from  a depth  x analyzed  at  the  angle 
9. 

Normalizing  the  intensity  with  respect  to  K,(0)  and  a. 


and  comparing  equation  (4-2)  with  the  Laplace  transform  of  a 
function  f(x), 


it  is  clear  that,  within  the  assumptions  of  the  model,  the 
normalized  XPS  intensity  as  a function  of  take-off  angle  is 
a Laplace  transform  of  the  composition  depth  profile  (CDP) 
of  the  specimen,  with  p = 1 / A sin(0)  as  the  parameter  of 
the  exponential  term. 

Several  studies  have  made  use  of  this  model  of  the 
generation  and  attenuation  of  the  photoelectron  flux  to 
extract  quantitative  information  concerning  inhomogeneities 
in  the  CDP.  For  the  most  part  the  samples  have  been  modeled 
as  a homogeneous  substrate  having  an  abrupt  interface  with  a 
single  homogeneous  overlayer.94'98 


(4-2) 


(4-3) 


These  single  overlayer 
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models  have  limited  applicability  in  surface  science.  It  has 
recently  been  demonstrated  that  reconstruction  of 
non-monotonic  CDPs  (i.e.  profiles  with  both  positive  and 
negative  gradients)  is  achievable  by  layer-wise  iterative 
solution  by  curve-fitting  using  Laplace  transforms.92  The 
fitting  curves  are  thus  self-consistent  with  the  exponential 
flux  attenuation  model.  However,  it  must  be  born  in  mind 
that  a unique  solution  with  good  quantitative  agreement 
between  depletion/accretion  depths,  magnitudes,  and 
gradients  cannot  be  uniquely  and  accurately  determined  from 
ARXPS  data.92 

The  model  flux  equation  assumes  an  atomically  smooth 
sample  surface,  and  negligible  dispersion  effects  on  both 
photon  and  photoelectron  fluxes.99  The  single  crystal 
substrates  analyzed  in  the  present  work  exhibit  minimal 
surface  roughness,  while  the  altered  surface  layers  are 
expected  to  have  an  amorphized  structure  that  will  wash  out 
the  diffraction  peaks  or  forward-scattering  effects  observed 
in  angle-resolved  photoemission  from  single  crystals.100 

Example  1:  Oxidized  Silicon 

The  Si  samples  were  analyzed  in  a Perkin-Elmer  PE5950 
XPS  spectrometer.  This  instrument  provides  high  angular 
resolution  and  high  count  rates  even  at  very  low  grazing 
angles.  The  resulting  uniform  angular  response  permits  data 
collection  across  a wide  range  of  sample  tilt  angles  [Figure 
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4-2].  The  spectrometer  was  operated  at  a resolution  of  0.01 
eV  full  width  at  half  maximum  (FWHM) , and  at  a pressure  of 
less  than  10'8  Torr.  The  Si  wafers  were  chemically  cleaned, 
then  oxidized  (in  the  ambient  atmosphere) . The  samples  were 
characterized  using  a 0 to  1000  eV  wide  scan  for  elemental 
analysis  followed  by  examination  of  the  Si  2p  lineshape. 

The  XPS  spectra  of  the  Si  2p  region  was  recorded  at 
takeoff  angles  of  10°,  30° , 60°,  and  90°.  After  background 
correction  the  peaks  were  curve-fitted  using  a PHI  software 
program.  In  addition  to  the  Si  2p  doublet  associated  with 
the  elemental  substrate,  shifted  peaks  associated  with  SiO, 
Si203,  and  Si02  bonding  were  resolved  [Figure  4-3,  Table 
4-1] . The  doublet  structure  was  not  resolved  in  the  shifted 
peaks.  The  resolved  peaks  were  integrated  and  normalized  to 
the  total  Si  2p  peak  at  each  angle,  resulting  in  a plot  of 
fractional  flux  versus  take-off  angle  for  each  suboxide 
[Figure  4-4].  These  data  were  then  solved  for  the  depth 
distribution  of  each  oxidation  state  independently,  without 
imposition  of  a mass  balance  constraint. 

The  solutions  [Figure  4-5]  show  the  Si02  as  the 
principle  constituent  of  the  oxide  overlayer  with  the  SiO 
and  the  Si203  at  the  Si02/Si  interface  in  the  expected 
sequence  of  decreasing  oxidation  with  depth.  The  composite 
of  these  separate  solutions  comes  close  to  unit  composition 
despite  the  lack  of  an  explicit  mass  balance  constraint.  The 
relative  intensities  of  the  all  three  surface  oxide  peaks  as 
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measured  at  the  normal  take-off  angle  fall  below  the 
expected  values,  while  the  elemental  Si  signal  is  stronger. 
This  is  understandable  as  the  effect  of  the  strong 
forward-scattering  peak  expected  from  the  substrate  along 
the  (100)  crystallographic  direction,  that  is,  the  surface 
normal . 96 

The  solutions  for  the  CDP's  are  plotted  as  functions  of 
relative  depth  x/A,  where  A is  the  mean  free  path  (MFP)  of 
the  Si  2p  photoelectrons.  The  value  of  A is  material- 
dependent;  it  can  be  taken  as  =27  ± 6 A in  the  Si  substrate, 
and  is  expected  to  increase  to  =37  ± 4 A in  the  oxide 
layer.96  Since  there  is  considerable  uncertainty  in  these 
values,  a constant  A is  used  in  these  plots  with  the 
understanding  that  the  effects  of  (1)  lower  number  density 
of  silicon  atoms  and  (2)  longer  electron  MFP  expected  in  the 
oxide  layer  will  act  in  opposition  in  the  flux  equation.  In 
essence  it  is  hoped  that  the  terms  cancel  out  in  plotting 
the  depth  profile  using  a constant  value  for  the  MFP. 

Example  2:  Sputter  Altered  Surfaces  of  (100)  GaAs 

The  GaAs  samples  were  analyzed  in  a Kratos  XSAM-800 
spectrometer.  The  wafers  were  cleaned  using  a 4:1:1  mixture 
of  H2S04:H202:H20.  They  were  then  Ar+  sputtered  using  a 
differentially  pumped  ion  gun  with  a rastering  beam  incident 
at  20°  relative  to  the  sample  surface.  During  the  sputtering 
the  system  pressure  was  held  at  10"9  Torr.  A steady-state 
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was  reached  at  each  of  the  ion  energies  (1.5,  3.0  and  5.0 
keV)  as  determined  by  monitoring  of  the  time  transients  in 
the  As  3d  and  the  Ga  3d  regions. 

ARXPS  spectra  from  the  GaAs  were  taken  at  takeoff 
angles  from  10°  to  90°  at  10°  intervals.  The  composition 
profiles  derived  [Figure  4-6]  show  a surface  composition 
near  the  bulk  value,  a subsurface  depletion  of  As,  then  a 
return  to  the  bulk  composition  at  greater  depth.  The 
subsurface  depleted  region  occurs  at  increasing  depth  with 
increasing  ion  energy.  The  MFP  of  the  sputter-altered  region 
has  been  assumed  equal  to  the  MFP  of  the  bulk. 

A review  by  Kelly101  described  a sputter-induced 
profile  similar  to  these,  with  a subsurface  depleted  region 
and  a surface  near  the  bulk  composition,  as  the  expected 
steady-state  for  an  alloy  exhibiting  bombardment-induced 
Gibbsian  segregation. 

Gibbsian  segregation  arises  from  a standard 
thermodynamic  treatment  of  a regular  solution  alloy,  having 
an  arbitrary  heat  of  mixing,  AHm,  but  an  ideal  entropy  of 
mixing.  If  there  is  an  enthalpy  change  associated  with  the 
exchange  of  an  A atom  on  the  surface  with  a B atom  in  the 
bulk,  then  the  steady-state  surface  composition  will  be101 


CA(s) 

CB(S) 


CA(b) 

CB(b) 


exp  (--§=,) 
RT 


(4-4) 


where  Q is  the  heat  of  segregation,  with  Cj(s)  as  the  surface 
concentration  and  C;,u,  as  the  bulk  concentration  of 
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component  i.  The  enthalpy  term,  Q,  may  be  attributable  to 
both  chemical  binding  and  strain  effects. 

Under  conditions  of  steady-state  sputtering,  the 
surface  composition  is  evaluated  from  the  conservation 
relation, 


Ca(s)  Ya  _ ^a(s)  y _ ^A(b) 

^B(s)  Yb  Cg(s)  C B(jb ) 


(4-5) 


where  Y(  is  the  sputter  yield  for  ions  impinging  on  the 
component  i in  the  target  and  YAB  is  a ratio  taking  into 
account  preferential  effects  such  as  energy  transfer  or 
chemical  binding.101 

In  Gibbsian  segregation  the  surface  composition  Cj(s)  is 
related  to  the  subsurface  composition  Cj(2)  according  to 
equation  (4-6) , that  is, 


CA(S) 

CB(s) 


CA(  2) 
CB(  2) 


exp  ( - -0- ) 
RT 


(4-6) 


Thus  during  sputtering  the  subsurface  layer  will  be  depleted 
of  the  component  which  segregates.  For  ion-bombarded 
specimens,  diffusion  in  the  undamaged  bulk  region  is 
expected  to  be  slow  compared  to  the  sputter  removal  rate.  A 
very  steep  concentration  gradient  will  result,  with  mass 
transport  largely  confined  to  the  ion  damage  region  where 
the  high  defect  density  will  lead  to  enhanced  diffusion 


rates. 
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This  mechanism  for  the  observed  preferential  sputtering 
and  the  resulting  composition  profile  help  to  explain  the 
conflicting  results  in  preferential  sputtering  of  different 
sys'tems  (e.g.  CuAu  alloys)  when  AES  vs.  ISS  has  been  used  to 
study  the  phenomenon. 102  Analysis  by  XPS  or  AES  with  the 
analyzer  at  a fixed  takeoff  angle  to  the  sample  surface  will 
indicate  a depletion  of  the  surface-segregating  species  due 
to  the  predominance  of  the  depleted  region  in  the  analyzed 
volume,  while  analysis  by  ISS  will  show  the  topmost  layer  to 
have  approximately  the  bulk  composition. 

The  GaAs  profiles  derived  from  ARXPS  support  this  model 
for  preferential  sputtering  effects.  The  increasing  extent 
of  the  depleted  region  with  increasing  ion  energy  suggests 
that  this  disturbed  layer  is  coincident  with  the  penetration 
range  of  the  bombarding  ions  for  these  energies. 

Conclusion 

ARXPS  has  been  demonstrated  as  a practical  technique 
for  non-destructive  depth  profiling  of  an  extremely  narrow 
surface  region  of  suitable  samples.  In  the  case  of  the 
oxidized  Si  samples,  the  location  and  relative  abundance  of 
Si  in  each  of  these  oxidation  states  was  recovered  from  the 
ARXPS  data.  The  sputtered  GaAs  composition  profiles  show  a 
surface  composition  near  the  bulk  ratio,  a subsurface  As 
depletion,  then  a return  to  the  bulk  composition  at  greater 
depth.  This  supports  the  mechanism  of  preferential 
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sputtering  due  to  Gibbsian  surface  segregation,  the  As  on 
the  surface  being  replenished  by  diffusion  from  the 
disturbed  layer  with  a diffusion  rate  enhanced  by  the  high 
defect  density  in  the  bombarded  region. 

The  application  of  these  results  to  this  study  gives  a 
more  detailed  understanding  of  the  surface  oxide  of  an  air- 
exposed  silicon  substrate,  for  example.  While  this 
contamination  layer  may  be  readily  sputtered  away  before 
deposition  of  the  silicide,  for  GaAs  sputter-cleaning  means 
the  silicide  will  be  deposited  on  a non-stoichiometric 
surface.  This  can  be  expected  to  have  an  influence  on  the 
interface  formation  and  contact  properties,  as  will  be 
discussed  in  Chapter  Six. 

Sputter-induced  changes  in  surface  composition  also 
adds  a complication  to  routine  surface  analysis  of  compounds 
such  as  the  silicide  films,  since  sputter  cleaning  and 
sputter-depth  profiling  are  commonly  used  in  Auger  electron 
spectroscopy  as  well  as  XPS.  The  alteration  of  surface 
composition  in  the  sample,  along  with  a number  of  other 
factors,  must  be  considered  in  a full  treatment  of  AES  data, 
to  be  taken  up  next  in  Chapter  Five. 


Table  4-1  Binding  energy  separation  [in  eV]  of  the  peaks 
deconvoluted  from  the  Si  2p  XPS  region  as  shown  in  Figure 
4-3.  The  elemental  2p  3/2  peak  is  used  as  a reference. 


Binding  energy  separation  [in  eV] 
measured  from  Si  2p  3/2  (elemental) 

Angle 

,Elem.  Si 

2p  1/2 

Si2+ 

Si3+ 

Si4+ 

90° 

0.67 

1.77 

2.67 

4.17 

30° 

0.67 

1.77 

2.67 

4 . 17 

20° 

0.67 

1.77 

2.47 

4.17 

10° 

0.63 

1.83 

2.73 

4.23 
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Surfac* 

Normal 


Figure  4-1  The  geometry  of  the  analyzer.  The  signal 
recorded  is  constrained  to  the  takeoff  angle  0,  and  the 
effective  path  to  the  surface  for  flux  from  a point  x below 
the  sample  surface  becomes  x/sin(0) . 
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SAMPLE  TILT  ANGLE  (deg) 


Figure  4-2  Plot  of  signal  intensity  versus  sample  tilt 
angle  for  the  Perkin  Elmer  PE5950  ESCA 
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BINDING  ENERGY,  «V 


Figure  4-3  Silicon  2p  spectra  from  a silicon  wafer  with  a 
native  oxide  layer,  at  three  sample  tilt  angles.  The  raw 
spectrum  has  been  Gaussian  curve-fitted  with  five  components 
to  the  fit:  the  Si  2p  elemental  peak  (doublet) , along  with 
three  shifted  peaks  associated  with  the  Si02,  SiO  and  Si203 
oxidation  states. 
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Figure 

versus 

angle. 


>.  Sine 


4-4  Plot  of  relative  silicon  species  intensities 
analysis  depth  Asin(0) , where  0 is  the  sample  tilt 
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Figure  4-5  ARXPS  results  for  the  oxidized  silicon  wafers, 
with  relative  intensity  measurements  [data  points  R(0)  = 

I (SiOx)/I  (SiT0T)  ] for  each  of  the  components  of  the  Si  2p 
region  of  the  XPS  spectrum,  plotted  versus  takeof f-angle  0 

(see  Figure  6-3) . Curve-fit  [ ] represents  the  Laplace 

transform  of  the  step-function  CDP  depicted.  The  atomic 
percent  plotted  is  actually  the  relative  percentage  of  Si  in 
each  oxidation  state,  ignoring  other  species  present  (such 
as  oxygen) . (a)  elemental  Si,  (b)  SiO,  (c)  Si203,  (d)  Si02, 

(e)  a composite  of  the  CDPs  for  all  subspecies.  The 
solution  profiles  are  consistent  with  a layered  structure  of 
sub-oxides  at  the  interface.  The  solutions  were  derived 
independently,  without  constraint  to  yield  100%  as  the  total 
composition. 
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Figure  4-6  ARXPS  results  for  the  Ar+  sputtered  GaAs , using 
the  Ga  3d  and  As  3d  XPS  peaks.  Relative  As  intensity 
(IAs/ (IAs+IGa) ) is  plotted  as  data  points  R(0)  , with  curve- 

fit  [ ] using  the  Laplace  transform  of  the  step-function 

CDP  depicted  (atomic  percent  As  vs.  relative  depth,  x/g) : 
for  Ar+  ion  energies  of  (a)  1.5  keV,  (b)  3.0  keV,  (c)  5.0 
keV.  Note  the  region  of  subsurface  As  depletion,  and  the 
progression  of  the  depleted  region  to  greater  depth  with 
increasing  ion  energy. 


CHAPTER  FIVE 

QUANTITATIVE  AUGER  ELECTRON  SPECTROSCOPY 
IN  THE  Gd-Si  SYSTEM 

Introduction 

Quantification  of  Auger  electron  spectra  from  rare 
earth  silicides  requires  correction  of  the  elemental 
sensitivity  factors  due  to  matrix  effects.  In  this  analysis 
the  Si  KLL,  Si  LW,  Gd  MNN,  and  Gd  NNN/NW  Auger  transitions 
have  been  measured  for  a number  of  samples.  Composition  was 
calculated  using  the  data  from  the  higher  kinetic  energy 
regions  (Si  KLL  and  Gd  MNN) . Comparison  of  Auger  data  with 
RBS  results  yielded  composition-dependent  effective  AES 
sensitivity  factors,  Seff,  using  the  RBS  stoichiometries  of 
the  codeposited  films  as  the  necessary  reference  bulk 
analysis. 

In  addition,  several  qualitative  or  semi-quantitative 
observations  are  made  including  (1)  rates  of  oxygen  uptake 
by  sputter-cleaned  surfaces  and  consequent  changes  in  the 
AES  lineshape,  (2)  changes  in  AES  lineshape  in  the  Si  KLL 
region  attributable  to  silicide  formation,  (3)  loss  features 
in  the  elastically  backscattered  primary  beam  and  comparison 
with  the  backscatter  from  elemental  samples  of  the 
constituents . 
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Experimental 

Samples  analyzed  included  thin  films  of  disilicide 
formed  by  reaction  of  Gd  metal  films  on  Si  substrates  and 
co— sputter- deposited  Gd— Si  mixtures.  In  order  to  confirm 
the  compositions,  a sample  of  bulk-prepared  disilicide  was 
polished  for  use  as  a standard  in  the  analyses.  Details 
concerning  the  RBS  and  AES  instrumentation  have  been  covered 
in  Chapter  Three. 

The  RBS  analysis  of  the  reacted  thin  films  was 
consistent  with  the  reported  single-phase  GdSi1  7 reaction 
product,  but  the  reaction  produces  a pitted  surface  and  a 
nonplanar  interface  with  the  substrate,  reducing  the  level 
of  confidence  in  the  stoichiometric  values  derived  from  the 
RBS.  The  disilicide  standard  on  the  other  hand  had 
inclusions  identified  as  Sic.  In  view  of  the  Sic  inclusions 
in  the  bulk  disilicide  and  the  roughness  of  the  reacted- 
thin-film  interface,  it  was  apparent  that  the  only  reliable 
calculated  RBS  stoichiometries  were  those  for  the 
codeposited  films. 

AES  composition  values  were  derived  from  the  data  using 
peak-to-peak  measurements  in  the  dN(E)/dE  mode  after 
numerical  differentiation.  The  dN(E)/dE  measurements  are 
strongly  affected  by  the  modulation  voltage  applied  to  the 
analyzer103  (or  eguivalently  the  differential  interval 
applied  to  the  N(E)  data).  In  this  work,  a differential 
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interval  equivalent  to  a 4 volt  peak-to-peak  modulation 
voltage  was  used  in  the  differentiation. 

The  RBS  compositions  of  the  codeposited  films  were 
found  to  fall  in  a range  of  Si:Gd  = 1.0  to  14.0.  Since  the 
composition  of  the  freely-formed  reaction  product  was 
regarded  as  desirable  in  terms  of  stability,  and  given  the 
difficulty  in  comparing  the  RBS  results  among  the  three 
types  of  samples,  it  was  necessary  to  develop  a capability 
to  assess  the  composition  of  the  codeposited  films  in 
comparison  with  regions  of  fully  reacted  thin  film  samples 
and  the  precipitate-free  matrix  of  the  bulk  disilicide. 

Thus  the  problem  was  defined  as  determining  the  AES 
fingerprint  of  the  GdSi1  7 disilicide  phase,  and  assessing 
the  effect  of  composition  and  contamination  on  AES 
measurements  in  the  Gd-Si  system. 

Quantitative  AES 


The  quantification  of  Auger  electron  data  begins  with 
an  expression  for  the  Auger  electron  current,  Ij(UVW) 
(hereafter  I . (Xf ) ) , from  the  UVW  Auger  transition  of  the  ith 
element104 


I?™=I  ■ 


UVW 


^(BP.ED)  c(Ep,En)  R 


■T-N(Xf)  'X(Xf)  - r[Ep,x’)  -Xf 


where  Ip=  primary  current,  T(UVW)  = Auger  transition 
probability,  <p.  (Ep,Eu)  = ionization  cross  section  of  the  U 
level  which  is  a function  of  the  primary  energy  Ep  and  the 
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critical  energy  for  U level  ionization  Eu,  Cj(Ep,Eu)  = the 
Auger  and  Coster-Kronig  contributions  to  electron  holes  in 
the  U level,  R = surface  roughness  factor,  T = instrument 
response  function,  N = atomic  density,  A = electron 
inelastic  mean  free  path  (IMFP) , r.  = electron  backscatter 
coefficient,  and  XjS  = atomic  fraction  of  the  ith  component 
in  the  volume  detected  at  the  surface. 

The  factors  T,  T,  , and  C,  are  matrix  insensitive.  The 
matrix  sensitive  factors  N,  A and  r are  functions  of 
concentration.  Ip  is  subject  to  instrumental  control,  within 
the  restrictions  on  its  accurate  measurement.  R may  vary 
from  sample  to  sample  depending  upon  initial  topography  and 
sputter-induced  changes  in  topography.  Lastly,  the  surface 
composition  XjS  may  differ  from  the  bulk  due  to  surface 
segregation  especially  under  the  ion  bombardment  that  is 
commonly  used  to  prepare  a clean  surface  for  analysis. 

Quantification  of  data  from  AES  can  be  attempted  using 
one  of  three  approaches: 105 

1.  First-principles  approach:  evaluation  of  terms  in 
the  Auger  current  equation. 

2.  Elemental  standards:  correcting  only  for  matrix 
effects . 

3.  Use  of  a series  of  standards  with  a composition 
near  that  of  the  unknown. 

The  first  two  have  the  great  advantage  in  that,  if 
successful,  they  implicitly  have  the  potential  for  allowing 
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samples  of  arbitrary  composition  to  be  quantitatively 
analyzed.  However,  in  actuality  the  evaluation  of  the  terms 
describing  both  matrix-insensitive  and  matrix-sensitive 
terms  is  dependent  upon  empirical  fitting  of  experimental 
measurements,  with  significant  scatter  in  the  fits  and  a 
resulting  large  margin  of  error  in  the  evaluations.105  On 
the  other  hand,  the  use  of  standards  minimizes  or  eliminates 
the  effects  of  many  factors  in  the  AES  flux  equation.  In 
addition,  by  comparison  of  results  from  AES  with  results 
from  a more  quantitative  technique  and  using  sensitivity 
values  from  elemental  standards,  one  can  derive  a single 
correction  term  S1eff  for  matrix  M that  will  bring  all  of  the 
results  into  agreement. 

This  can  be  expressed  as  follows:  for  a pure  elemental 
sample,  an  elemental  sensitivity  factor  Sj°  can  be 
defined106-107  as 


where  Ij0  is  the  Auger  current  from  a pure  sample  of  element 
i and  Ip  is  the  primary  beam  current. 

For  a pure  elemental  sample  there  can  be  no  segregation 
effects,  so  Xsi  = X. . Therefore  from  Eqs.  (5-1)  and  (5-2), 

S/*rTOr-V/'-V  ■ciE'-E",-R-T-Ni°-X°rri°(E„)  . 


(5-3) 
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From  this  expression  can  be  derived  the  relative 
elemental  sensitivities  Sjrel  as  tabulated  in  the  Handbook  of 


where  the  subscript  S refers  to  a standard  reference  peak  in 
a pure  elemental  sample  (typically  the  Ag  351  eV  peak) . The 
matrix-insensitive  factors  relating  to  ionization  and 
transition  probabilities  have  been  condensed  into  the  terms 
G.  and  Gs,  and  differences  in  surface  roughness  are  ignored. 

Presuming  that  the  natural  shape  of  a given  feature  is 
described  by  G( , and  that  the  analyzer  and  detector  response 
during  the  acquisition  (including  modulation  voltage  or 
numerical  differentiation  window  for  dN ( E ) / dE  measurements) 
is  the  same  as  that  for  the  elements,  correct  measurements 
then  depend  upon  evaluation  of  the  terms  involving  the 
matrix-dependent  factors  N,  A and  r.  Thus  all  the  elemental 
parameters  are  normalized  in  a self-consistent  fashion,  with 
the  caution  that  the  transmission  of  the  analyzer  and 
detector  response  may  vary  greatly  under  various  operating 
conditions  of  the  same  instrument. 


Auger  Electron  Spectroscopy108 


(5-4b) 


( 5-4a) 


(5-4c) 
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Matrix-Corrected  Sensitivity  Factor  Seff 


Relative  elemental  sensitivities  can  in  many  cases  be 
evaluated  directly  by  comparison  of  measured  intensities 
from  prepared  standards  (Eq.(5-4b)),  without  evaluating  the 
separate  terms  of  Eq. (5-4c) . However,  for  evaluating  the 
composition  of  arbitrary  samples  a more  general  approach  is 
desired. 

To  parallel  the  notation  used  for  the  matrix  terms  by 
Mathieu  and  Landolt,109  a corrected  sensitivity  factor  Sieff 
for  the  elemental  ransition  i in  matrix  M can  be  expressed 
as  the  product  of  the  elemental  sensitivity  S,.0  and  a matrix 
correction  factor  k^- : 


(5-5) 


where 


(5-6a) 


and 


( 5-6b) 


(5-6e) 


(5-6d) 


(5-6C) 


( 5— 6f ) 
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The  superscript  M refers  to  values  pertaining  to  the  given 
matrix;  the  superscript  0 refers  to  the  elemental  matrix  i. 

Combining  Eqs. (5-1) -(5-3)  and  Eq. (5-5)  gives  the 
intensity  from  element  i in  matrix  M, 


For  comparison  of  peaks  i and  j of  two  elements  in  the 
matrix  M, 


Evaluation  of  matrix-dependent  terms 

Attempts  have  been  made110'115  to  evaluate  the  matrix- 
dependent  terms  in  Eq. (1)  which  bear  directly  on  the  terms 
k,M,  kjA,  and  kjr,  describing  the  number  density  of  atoms,  the 
electron  mean  free  path,  and  backscatter  of  the  primary 
electron  beam.  Hall  and  Morabito110  evaluated  the  matrix- 
dependent  parameters  in  detail,  with  a focus  on  dilute 
alloys. 

(The  terms  corresponding  to  segregation  and  analyzer 
transmission,  kjS  and  kj,  are  not  dealt  with  in  these 


(5-7) 


and  the  atomic  fraction  of  element  i 


(5-8a) 


( 5-8b) 


82 


treatments,  but  as  noted  they  can  have  a serious  impact  upon 
results) . 

For  the  backscatter  term,  Hall  and  Morabito110  relied 
upon  an  empirical  equation  given  by  Reuter,111  designed  for 
use  with  the  electron  microprobe: 

r/=  1 + 2.8  (1-0.9  Ec/Ep)  Q7,  (5-10) 

where  Ec  is  the  critical  ionization  energy  of  the  initial 
core  hole  and  O'  is  the  matrix-dependent  "backscatter 
coefficient",  given  for  pure  elements  by 

Q/  = -0.0254  + 0.016Z  - 0 . 000186  Z2  + 8.3xlO_7Z3  (5-11) 

where  Z is  the  atomic  number  of  the  matrix  element. 

On  the  other  hand,  Prutton  and  El  Gomati112  have 
suggested  that  the  backscatter  factor  is  linear  in  atomic 
number,  based  upon  analysis  of  spectral  backgrounds.  The 
straight  line  fit  given  is 

r = 1 + 0.15  +0.152Z  . (5-12) 

For  the  inelastic  mean  free  path  A,  Hall  and  Morabito110 
considered  both  the  results  of  Penn113  and  those  of  Seah  and 
Dench.114  While  for  a given  element  the  results  derived  from 
these  two  approaches  may  be  in  considerable  disagreement,  in 
general  the  values  are  similar. 

There  is  an  advantage  in  using  the  formula  for  A from 
Seah  and  Dench, 

A = A-N-Ui-E~2  + B(E/N)1/2,  (5-13) 

which  for  higher  energies  reduces  to  the  second  term.  Then 
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with  Eq.  ( 5-6c)  above  k,A  can  be  approximated  as 

ki  = (Nf/N?)1'2  . 


(5-14) 


Number  density  in  an  unknown  sample  may  be  estimated, 
or  else  approximated  by  weighting  the  elemental  densities 
according  to  fractional  composition. 

In  this  work,  a modified  version  of  the  Reuter  term  was 
used  in  which  calculated  elemental  backscatter  values  were 
weighted  according  to  the  sample's  fractional  composition. 
Comparable  results  were  obtained  using  a simpler  expression 
from  Prutton  and  El  Gomati,112  which  will  be  used  below. 

The  term  for  number  density  was  estimated  by  a weighted 
average  of  the  atomic  number  density  of  the  silicide  and  the 
elemental  density  of  the  element  that  is  in  excess,  assuming 
formation  of  silicide  is  maximized.  These  values  were  found 
to  be  consistent  with  measured  values  for  density  as 
determined  by  calculations  based  on  profilometry 
measurements  of  film  thickness  combined  with  measured  mass 
gain  on  a crystal  thickness  monitor. 

Combining  (1)  an  evaluation  of  kt.N,  (2)  an  expression 
for  k/,  also  in  terms  of  number  density  (Eq. (5-14)),  and 
(3)  evaluating  k.r  using  Eq.  (5-12)  with  ZM  as  the  average  Z 
for  the  matrix,  gives  an  evaluation  of  the  matrix  correction 
factor 


/#  (1.15  + 0.0152  ZM) 
(1.15  + 0.0152  Z° 


<«„/  w/> 1/2 


(5-15) 
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where  k'  is  introduced  as  a term  to  accommodate  errors  due 
to  the  approximations  and  other  factors,  such  as  possible 
surface  segregation  effects.  It  is  emphasized  once  again 
that  the  factor  for  analyzer  transmission,  kj  (Eqs. (5-6a) , 

( 5-6f ) ) can  lead  to  significant  variations  in  measured 
intensity,  given  the  warnings  by  Powell115  et  al.  concerning 
variations  in  Auger  spectra  between  laboratories.  It  is 
apparent  that  variations  in  the  operating  conditions  even  on 
the  same  instrument  may  have  significant  impact  on  results 
as  well.  In  the  present  work  this  was  dealt  with  explicitly 
by  taking  measurements  from  elemental  standards  including 
comparison  with  the  Ag  351  eV  peak  before  applying  the  above 
correction  term,  and  thus  ksT  is  treated  as  unity. 

The  value  of  this  approach  is  twofold.  It  allows  the 
immediate  development  and  employment  of  quantitative  AES 
analysis  of  samples  within  the  range  of  certainty  of  the 
standards,  with  results  that  represent  the  greatest  degree 
of  accuracy  currently  obtainable  from  AES.  The  advantages  of 
AES  in  rapid  data  acquisition,  good  detection  limits, 
relatively  uniform  elemental  response,  and  lateral 
resolution  combined  with  reasonably  accurate  quantification 
can  then  be  employed  to  great  advantage  in  an  ongoing  study, 
without  the  need  for  repetitive  compositional  analyses  by 
the  more  expensive  and  time-consuming  quantitative 
techniques . 
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Secondly,  after  applying  the  correction  term  as  a 
function  of  composition  it  is  reasonable  to  discuss  which 
factors  in  the  Auger  current  equation  would  serve  as 
reasonable  explanations  for  any  residual  disagreement 
between  AES  and  bulk  analysis  (as  expressed  by  k'  in  Eq. (5- 
15)).  While  reasonable  values  for  number  densities  may  be 
available,  there  remains  considerable  uncertainty  in  the 
evaluations  of  A and  the  backscatter  r.  In  particular,  the 
assessment  of  r ignores  the  differences  to  be  expected  in 
the  susceptibility  of  peaks  having  different  initial  hole 
ionization  energies,  with  respect  to  a backscatter 
excitation  flux  that  would  be  expected  to  vary  considerably 
as  a function  of  energy.  Assessments  of  this  nature  will 
ultimately  serve  to  increase  our  knowledge  of  the  first 
principles  involved,  and  thus  bring  us  closer  to  the 
realization  of  quantifiable  AES  of  arbitrary  unknown  samples 
using  the  first-principles  approach. 

In  summary,  the  strategy  is  now  defined  as  evaluation 
of  composition  using  simple  corrections  to  elemental 
sensitivity  factors  which  have  been  corrected  for  the 
response  of  the  instrument.  Comparison  of  these  results  with 
RBS  results  give  values  of  k'  that  allow  us  to  assess  the 
accuracy  of  the  model-derived  parameters  for  N,  A and  r.  It 
is  then  feasible  to  discuss  other  possible  factors  in  the 
AES  data  analysis,  for  instance  sputter-induced  surface 
composition  changes  as  reported  by  Wirth  et  al.116  in 
quantitative  AES  analysis  of  refractory  metal  silicides. 
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Results  and  Discussion 


Qualitative  Results 

Elastic  peak  and  backscatter 

The  elastically  backscattered  primary  beam  and  its 
background  was  studied  for  several  elemental  samples  and  the 
silicides  by  lowering  the  energy  Ep  of  the  primary  beam  to 
2.8  keV,  to  fall  within  the  limits  of  the  detector  [Figure 
5-1 (a) ] . The  intensity  of  the  elastically-scattered  primary 
beam  along  with  the  background  level  at  lower  energies  in 
these  spectra  provide  a qualitative  measure  of  the 
backscatter  factor,  and  also  indicate  characteristic  loss 
peaks  for  each  material.  The  oxidized  surface  of  the 
silicide  exhibits  a strong  loss  feature  at  =30  eV  below  the 
elastic  peak,  and  significant  reduction  of  the  elastic  peak 
height  compared  to  the  clean  silicide  surface.  As  the  oxide 
was  removed  by  ion  sputtering,  the  loss  peak  at  30  eV 
decayed,  while  two  new  loss  peaks  appeared  at  =12  eV  and  =27 
eV  [Figure  5-1 (b) ] . 

Low  energy  region 

In  the  interval  of  0-160  eV  occur  the  Si  LW  and  Gd 
NNN/NW  Auger  lines.  Because  of  the  higher  surface 
sensitivity  due  to  the  short  IMFP  for  lower  energy  electrons 
and  a nonlinear  response  of  the  detection  system  in  this 
energy  regime,  these  peaks  are  not  optimal  for  purposes  of 
quantification.  The  well-known  chemical  shift  in  the  Si  LW 
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peak  associated  with  oxide  formation  is  also  a negative 
factor  in  this  regard  [Figure  5-2],  but  the  observed  change 
in  lineshape  does  serve  as  a useful  monitor  of  oxygen 
incorporation  in  the  silicides  and  oxygen  uptake  from  the 
residual  gas  in  the  chamber.  Changes  in  the  112-127  eV 
portion  of  the  Gd  lineshape  were  also  observed  on  the 
oxidized  surface,  with  the  112  eV  and  the  127  eV  peaks 
decaying  along  with  the  appearance  of  a peak  at  122  eV.  The 
138  eV  Gd  NW  peak  is  not  noticeably  affected  by  this 
oxidation  [see  Figure  5-2]. 

High  energy  regions 

Two  regions  of  higher  energy  peaks  were  studied.  The 
Gd  MNN  peaks  occur  between  600-1250  eV,  while  the  Si  KLL 
peaks  are  found  at  1520-1640  eV.  By  observation  of 
lineshape  changes  these  peaks  were  found  to  be  not  strongly 
affected  by  adventitious  oxygen  in  the  chamber.  It  was  in 
these  regions  that  reliable  quantitative  results  were 
obtained. 

A pronounced  change  in  the  Si<LL  lineshape  was  observed 
upon  formation  of  the  rare  earth  silicide  [Figure  5-3].  The 
1604  eV  peak  in  that  region  decays  monotonically  with 
increasing  Gd  content.  This  feature  is  believed  attributable 
to  a plasmon  loss  mechanism,  and  the  changes  in  the 
lineshape  therefore  are  ascribed  to  changes  in  that  loss 
mechanism  due  to  incorporation  of  Gd. 
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Quantitative  Results 

Evaluation  of  k1  (residual  error) 

After  applying  the  calculated  corrections  to  the 
dN(E)/dE  peak-to-peak  intensities,  there  remained 
significant  disagreement  between  the  AES  results  and  the  RBS 
[Figure  5-4].  Angle-resolved  X-ray  photoelectron 
spectroscopy  (ARXPS)  data  indicated  that  sputter  cleaning 
creates  a silicon-rich  surface  layer  of  about  loA.  However, 
the  effect  on  the  AES  data  should  be  minimized  since  the 
flux  entering  the  wide  acceptance  angle  of  the  AES  analyzer 
is  dominated  by  the  signal  from  near  the  surface  normal. 
Besides,  direct  comparison  with  AES  data  from  samples  of  the 
elements  shows  not  only  silicon  but  also  Gd  at  higher 
concentration  than  expected. 

Normalizing  the  independent  results  to  100%  total 
composition  (i.e.,  eliminating  the  factor  of  beam  current 
from  Eq.  5-1)  does  not  improve  the  result  due  to  the 
countervailing  senses  of  the  deviations  for  the  Si  and  the 
Gd  intensities  (k'sj,  k'Gd)  . Thus  it  is  reasonable  to  discuss 
this  disagreement  between  the  RBS  and  the  AES  in  terms  of 
the  model  factors  N,  A and  r. 

Linear  fit  to  experimental  deviation  (k11) 

The  residual  error  k'  was  plotted  against  the  corrected 
AES  atomic  percent,  and  empirically  fitted  with  a linear 
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equation  [Figure  5-5].  For  the  Gd  data,  the  residual  error 
was  compensated  by  a factor  k"  given  by 

*&-  2.M.KC"']  <5-16> 

and  for  Si, 

JCsi  = 0.2  +0.8  (x£S'corz)  . (5-17) 

While  a more  detailed  fit  to  the  data  is  feasible,  the 
large  deviation  of  the  error  factor  from  unity  (k'  varies 
from  0.4  to  2.8)  means  that  even  this  very  simple  approach 
to  a correction  term  will  greatly  improve  the  results. 
Results 

The  results  were  encouraging  [Tables  5-1  and  5-2, 

Figure  5-4]  for  films  ranging  from  50%  Si  to  100%  Si,  but 
the  samples  near  the  1.7:1  stoichiometry  of  the  disilicide 
phase  were  not  brought  into  agreement.  This  group  of  samples 
included  the  bulk  disilicide,  the  reacted  thin  films,  and 
the  codeposited  disilicide.  It  was  noted  that  the  deviation 
of  these  samples  were  indicative  of  excessive  Gd  and  also 
excessive  Si  [see  Figures  5-4  and  5-5] ; though  the 
independent  assessment  of  atomic  percent  of  each  element 
gave  a total  greater  than  100%,  the  ratio  of  Si:Gd  was  very 
near  to  the  expected  1.7:1. 

When  the  percentages  were  normalized  to  total  100% 
there  was  good  agreement  with  the  RBS  analysis  throughout 
the  range  of  compositions,  including  the  anomalous 
disilicides,  with  final  results  in  agreement  within  1.5 


Atom%  RMS  error. 
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Discussion 

The  ratios  of  the  low-energy  peaks  to  the  higher-energy 
peaks  (i.e.,  Si  LW/  Si  KLL  and  Gd  NNN/Gd  MNN  [Figure  5-6] 
provide  evidence  that  the  mean  free  paths  of  the  Auger 
electrons  is  strongly  affected  by  the  composition  of  the 
sample.  Incorporation  of  gadolinium  reduces  the  intensity 
of  the  silicon  peaks,  with  a more  severe  attenuation  of  the 
high  energy  Si  KLL  peak  compared  to  the  Si  LW.  The 
intensity  of  the  Gd  NNN  peak  is  augmented  to  a greater 
degree  than  the  Gd  MNN. 

These  observations  may  be  understood  as  resulting  from 
a shift  in  the  electron  mean  free  path  versus  kinetic  energy 
curve.  For  Si,  the  inclusion  of  Gd  with  its  numerous 
shallow  core  electrons  into  the  matrix  produces  an  abundance 
of  scattering  opportunities  with  respect  to  elemental 
silicon,  with  proportionately  more  of  these  energy-loss 
interactions  becoming  available  for  the  higher  energy  KLL 
electrons.  For  Gd,  the  substitution  of  Si  atoms  reduces  the 
possibilities  of  inelastic  scattering,  and  similarly  the 
high-energy  Gd  MNN  peaks  have  more  potential  loss- 
interactions  removed  by  the  substitution. 

The  anomalous  AES  response  of  the  disilicides  could 
have  arisen  from  a number  of  factors  such  an  increase  in  the 
primary  beam,  which  was  measured  by  biasing  the  samples.  But 
the  close  grouping  of  samples  prepared  using  very  different 
techniques  would  seem  to  indicate  an  intrinsic  property  of 
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the  disilicide.  The  AES  response  for  the  disilicides  is 
consistent  with  an  increase  of  =20%  in  the  mean  free  path  of 
the  electrons,  due  to  formation  of  a well-ordered  lattice  in 
a narrow  composition  range  with  consequent  reduced 
scattering  of  the  photoelectrons. 

Conclusion 

The  introduction  of  AES  equipment  with  good  beam 
stability  and  areal  resolution  along  with  digital 
acquisition  of  E'N(E)  spectra  provides  an  opportunity  to 
develop  reliable  quantitative  AES  analysis  in  simple 
systems.  Using  standard  samples  of  both  elements  and 
compounds  (whose  compositions  are  determined  by  some  method 
of  bulk  analysis) , matrix-sensitive  effective  sensitivity 
factors  have  been  derived  that  provide  the  means  to  correct 
the  Auger  current  measurements  and  therefore  determine 
compositions  quickly  and  accurately.  Qualitative  and  semi- 
quantitative  information  on  chemical  bonding  is  also 
realized  from  changes  in  peak  shape  in  the  N(E)  mode. 

Results  in  the  Gd-Si  system  demonstrate  that 
empirically  derived  linear  relationships  based  on  atom%  as 
measured  by  AES  can  yield  quite  accurate  effective 
sensitivity  factors.  It  is  suggested  that  the  measured 
variation  of  the  effective  sensitivity  may  be  readily 
attributable  to  one  of  the  first-principle  terms  in  the 
model  equation,  and  that  compilations  of  such  results  may 
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lead  to  better  understanding  of  AES  from  the  first- 
principles  approach. 

In  the  case  of  the  Gd-Si  system,  the  departure  from 
elemental  sensitivities  is  such  that  it  may  be  most  readily 
attributable  to  changes  in  the  inelastic  mean  free  path  A, 
with  differing  effects  for  the  Auger  energies  of  the  several 
transitions  investigated.  The  singularity  observed  in  the 
AES  response  of  the  disilicides  is  consistent  with  a marked 
increase  in  the  mean  free  path  of  the  photoelectrons  in  the 
single-phase  disilicide  compared  to  mixtures.  Further  work 
will  be  required  to  substantiate  this.  However,  in  the  range 
of  composition  of  this  study,  quantified  compositional 
information  is  now  readily  available  from  AES. 


Table  5-1  Atomic  percent  silicon  and  matrix  correction  factors. 

Atomic  percent  calculated  from  AES  data  (Si  KLL  peak) , both  the  raw  peak-to-peak 
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KINETIC  ENERGY  (eV) 


Figure  5-1  Elastic  backscattered  beam  (Ep  = 2.8  keV)  . 

(a)  Elastic  peak  for  various  samples  (Gd,  Ag,  GdSi1  7,  GdOx, 
Si) ; (b)  Transients  in  the  elastic  peak  and  characteristic 

losses,  during  sputter  cleaning  (3  keV  Ar+)  of  air-exposed 
GdSi1  7.  Note  the  immediate  increase  in  the  elastic  peak 
intensity  upon  removal  of  the  oxygen  contamination  and  the 
development  of  characteristic  loss  peaks  at  =12  eV  and  =27 
eV  below  the  primary  peak. 
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Figure  5-2  Effects  of  oxygen  contamination  on  gadolinium 
silicide  films  [Si  LW/Gd  NNN  AES  region] 

(a)  (1)  clean  surface  (2)  after  long-term  exposure  to 
residual  gas  in  UHV  chamber. 

(b)  (1)  clean  film  (2)  film  high  in  incorporated  oxygen 
during  deposition  (=20%) 
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KINETIC  ENERGY  (eV) 


At%Si  (RBS) 


Figure  5-3  Si  KLL  AES  region:  (a)  Comparison  of  spectra 
from  elemental  Si  and  from  GdSi1  7 show  a marked  quenching  of 
satellite  peak  'a',  (b)  Fractional  intensities  of  minor 
peaks  in  the  Si  KLL  region.  Peak  'c'  is  the  KL1 L1 
transition;  'a',  'b'  and  ' d'  are  loss  peaks.  Note  the 

monotonic  decay  of  the  loss  peak  'a'  with  increasing  Gd 
content. 
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PERCENT  Si  [RBS] 


Figure  5-4  AES  composition  of  GdSix  versus  RBS 
measurements:  (AES  raw  peak-to-peak,  matrix  corrected, 
corrected  by  empirically  derived  term  k'1). 

(a)  Si  KLL  (b)  Gd  MNN 


and 
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Figure  5-5  Error  factor  after  application  of  matrix 
correction  terms,  and  empirical  fit  k'1: 

(a)  Si  KLL  (b)  Gd  MNN 
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PERCENT  Si  [RBS] 


Figure  5-6  Comparison  of  high  energy  and  low  energy  AES 
intensities:  (a)  Si  KLL  and  LW  (b)  Gd  MNN  and  NNN.  With 
increasing  Gd  content,  the  lower  energy  transitions  are 
proportionately  higher  in  intensity  compared  to  the  higher 
energy  peaks. 


CHAPTER  SIX 
RESULTS  OF  ANALYSIS 


Crystallization  Behavior  of  Codeposited  Filins 
X-rav  Diffraction 

Analysis  using  Cu  KCT  radiation  indicated  that  the  films 
have  an  amorphous  character  as  deposited  and  crystallize  to 
a disilicide  structure  after  a brief  anneal  at  T=400°C.  The 
diffraction  pattern  then  obtained  showed  a better  agreement 
with  the  indexed  pattern117  for  orthorhombic  "GdSi2"  than  did 
the  reacted  thin  films  [Figure  6-1] , indicating  a more 
homogeneous  phase  formation  in  the  codeposited  thin  films. 

A strong  reflection  at  20=34°  (d=2.63A)  assigned  to 
orthorhombic  GdSi2  (112)  indicated  a high  degree  of 
texturing  in  the  polycrystalline  silicide. 

Thicker  samples  (5,OOoA)  from  a well-seasoned  sputter 
target  were  deposited  to  increase  the  diffracted  signal 
intensity.  These  produced  very  interesting  results  in  XRD, 
showing  only  two  very  dominant  features  even  before 
annealing  [Figure  6-2,  Table  6-1],  one  near  the  familiar 
GdSi2  (112)  peak  at  34°,  and  another  peak  not  previously 
observed  near  20=65°  (d=1.44A)  which  could  be  assigned  to 
GdSi2  (220) . The  dominance  of  these  peaks  is  indicative  of 
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two  distinct  strong  preferred  orientations  in  the  films. 
Similar  results  were  obtained  with  both  Si (100)  and 
GaAs(lOO)  substrates. 

Rapid  annealing  (500°C/20  sec.)  caused  these  two  peaks 
to  shift  to  lower  angles  (larger  d-spacings,  1.45A  and  2.65A 
respectively) . This  was  consistent  with  a lattice  expansion 
due  to  the  excess  silicon  in  the  deposition  (Gd/Si=2) . In 
fact,  the  ThSi2-type  tetragonal  phase  of  GdSi283  would  have 
spacings  of  2.67A  and  1.45A  for  the  (112)  and  (220) 
reflections  [Table  6-2],  and  it  may  be  nucleated  as  a 
metastable  form  under  these  conditions. 

Annealing  at  progressively  higher  temperatures  caused 
these  peaks  to  shift  back  toward  larger  20  values,  with  d- 
spacings  reduced  to  values  below  the  as-deposited  values. 

At  temperatures  above  700° C these  spacings  are  close  to  the 
values  calculated  for  the  (lTl)  and  (112)  planes  in 
hexagonal  A1B2  GdSi,  7 [see  Table  6-2].  Considering  the 
crystallographic  relationship  between  the  orthogonal  or 
tetragonal  half  unit  cell  and  the  hex  unit  cell 
[see  Figure  2-1,  above],  a phase  change  to  hexagonal 
structure  would  exhibit  related  diffraction  peaks  from  these 
closely  related  hexagonal  planes  [Figure  6-3]. 

The  hexagonal  GdSi1  7 phase  as  a significant  constituent 
of  thin  film  specimens  has  previously  been  observed  only  for 
liquid  phase  reaction  of  Gd  on  Si (111)  under  rapid  heating 
by  electron  beam.11,12  Thus  its  detection  in  these  films  on 
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Si (001)  is  surprising,  and  should  be  understood  as  tentative 
without  other  corroborating  diffraction  lines.  The 
similarities  in  the  structures  and  the  variation  in  reported 
lattice  parameters  from  bulk-sample  studies  adds  to  the 
uncertainty  [see  above,  Table  2-6]. 

The  relative  heights  of  the  two  peaks  varies  with  time 
and  temperature  [Figures  6-4  to  6-7].  For  rapid  anneals  (20 
sec.  RTA)  the  65°  peak  is  quenched  dramatically  compared  to 
the  as-deposited  film,  then  grows  slightly  with  temperature, 
doubling  from  600°C  to  800°C  but  remaining  smaller  than  as- 
deposited  [Figure  6-4,  (a)-(e)].  With  longer  times  (40 

sec.)  the  intensity  of  this  peak  increases  markedly  [Figure 
6-4,  (f)  and  (g) ; Figure  6-5]  and  nears  the  intensity  of  the 
34°  peak. 

The  34°  peak  grows  steadily  with  time  and  with 
temperature  [Figure  6-6]  and  begins  to  develop  a slight 
shoulder  at  33.3°  (d=2.69A)  that  becomes  a prominent 
satellite  at  800°C.  This  satellite  is  even  more  prominent 
in  the  as-deposited  film  but  is  totally  quenched  by  the 
500° C treatment. 

Furnace  treatments  (375°,  30  min.)  [Figure  6-7]  of 
films  on  GaAs  substrates  result  in  a single  major  peak  at 
34°  assigned  to  orthorhombic  GdSi2(112),  and  only  a small 
shoulder  on  the  substrate  peak  near  64°.  On  silicon  the  two 
peaks  occur  at  larger  angles  and  show  nearly  equal 
intensity,  more  similar  to  the  longer-term  RTA  heating  [see 
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Figure  6-5 (b) ] . After  1 hr.  at  900°C  [Figure  6-7 (d) ] the 
33.3°  satellite  peak  seen  in  the  800°C  RTA  data  dominates  a 
broad  amorphous-type  pattern.  No  clear  assignment  could  be 
made  for  this  d-spacing. 

Treatment  at  higher  temperatures  (near  800°C)  [Figure 
6-4]  and  for  longer  times  [Figure  6-4,  6-5]  also  produced 
broad  features  at  =51°,  d=1.79A,  which  may  be  assigned  to 
orthorhombic  GdSi2  (211)  and  (121),  d=1.78A  and  d=1.8oA 
respectively. 

In  summary,  the  strong  orientation  of  these  films, 
especially  the  single  strong  (112)  orientation  observed  on 
the  GaAs , indicate  that  the  silicide  may  form  a highly 
ordered  interface  with  semiconductor  substrates.  At  the 
same  time  the  phase  changes  implied  by  these  data  and 
supported  by  the  complex  crystallography  of  the  silicides 
show  that  much  remains  to  be  done  before  a clear 
understanding  will  be  reached. 

Transmission  Electron  Microscopy 

Attempts  to  clarify  the  XRD  data  using  electron 
diffraction  were  not  very  fruitful.  Samples  were  prepared 
for  transmission  electron  microscopy  (TEM)  in  both  the 
planar  and  cross-sectional  view,  revealing  a columnar  grain 
structure  with  grain  diameters  of  =50oA  [Figures  6-8  to  6- 
10],  Electron  diffraction  produced  a very  spotty  ring 
pattern  (polycrystalline) , with  several  bright  spots 
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associated  with  several  individual  grains  of  closely  related 
orientation. 

The  cross-sectional  samples  produced  had  very  few 
transparent  areas  due  to  difficulty  in  thinning  the  silicide 
without  destroying  the  substrate  support.  The  high  mass 
absorption  of  the  gadolinium  also  requires  the  sections  to 
be  quite  thin.  The  few  areas  available  were  imaged  along 
the  Si (Oil)  and  (111)  directions  in  hope  of  finding  silicide 
crystals  with  related  orientation  and  planes  parallel  to  the 
Si (001)  surface  [Figure  6-8,  6-10]. 

The  columnar  grains  do  not  seem  to  be  as  singularly 
oriented  as  indicated  by  the  XRD  data.  Several  large  grains 
were  correlated  with  plane  spacings  consistent  with  the 
lattice  spacing  of  hexagonal  GdSi2(lTl),  d=2.6lA  [Figure  6- 
9(d),  6-10 (b) ] . Some  of  these  grains  also  showed  strong 
contrast  from  other  distinct  planes  in  nearly  the  same 
direction,  with  a spacing  of  2.00A  [Figure  6-10 (c) ] . This 
result  is  not  well  understood,  but  could  be  due  to  the 
interpenetration  of  silicon  and  gadolinium  sublattices  with 
different  scattering  power.  Other  grains  with  lattice 
spacing  3.37A  could  be  assigned  to  hexagonal  GdSi2(110)  , 
d=3 . 37A  [Figure  6-9 (c) ] . 

The  discrepancy  between  these  data  and  the  XRD  results 
can  be  tentatively  explained  by  supposing  the  XRD  data  are 
generated  by  a small  number  of  perfectly  oriented  grains. 
With  the  rest  of  the  columnar  grains  only  slightly  out  of 
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alignment  with  the  silicon,  very  few  other  reflections  would 
be  detected  by  XRD.  The  data  are  thus  dominated  by  the 
oriented  fraction,  with  the  nearly  oriented  grains 
essentially  transparent  to  the  analysis. 

While  this  explanation  is  largely  conjectural,  no  other 
answer  offers  itself  and  further  investigation  will  be 
necessary  to  resolve  the  question. 

Metallurgical  Stability 

Stabilitv/Reactivitv  of  Films  on  Silicon 

Dramatic  evidence  of  the  role  of  silicon  as  the 
diffusing  species  in  the  rare  earth/silicon  reaction60  was 
observed  in  the  codeposited  films.  In  silicon-rich 
depositions,  moderate  anneals  (RTA  400-600°C,  20  sec.)  led 
to  extensive  recrystallization  of  the  GdSi1  7 phase,  with  the 
excess  silicon  diffusing  through  a 5,000A  film  to  the 
surface  where  it  aggregates  [Figure  6-11]. 

Furnace  anneals  up  to  900°C  for  one  hour  produced  a 
massive  metallurgical  reaction  [Figure  6-12],  with  silicon 
whiskers  growing  out  of  ruptures  in  a blistering  surface 
oxide.  This  presents  an  upper  limit  to  the  temperature 
stability  that  has  not  previously  been  observed.  The 
GdSi1  7/Si  interface  has  been  reported  to  be  stable  up  to 
900°C.9 

The  low  temperature  of  activation  for  the  rare  earth 
silicide  reaction  on  silicon45,57  explains  the  continued 
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reactivity.  It  appears  that  silicide  remains  intact  as  a 
catalytic  interface  while  silicon  continues  to  diffuse 
through  the  film  and  nucleate  and  grow  at  the  surface  in  a 
massive  diffusion.  The  relatively  high  heat  of  formation  of 
the  silicide4'118  (comparable  to  USi2,  AHf=  -31.2  kCal/mole) 
[Table  6-3]  would  prevent  its  dissolution,  while  the 
catalytic  effect  of  the  rare  earth  metal  atoms  on  silicon 
bonds  in  the  substrate  continues  to  provide  free  silicon. 

The  defect  lattice  of  the  silicide  provides  a large  number 
of  silicon  vacancies  to  enhance  diffusion  rates  through  a 
what  is  already  a relatively  open  lattice  structure. 
Fortunately,  these  conditions  are  beyond  device  performance 
limits,  but  they  are  conditions  that  could  be  approached  in 
processing  steps  and  further  investigations  are  called  for. 
Stabilitv/Reactivitv  of  Films  on  GaAs 

The  recrystallization  of  silicon-rich  films  on  GaAs  is 
similar  to  the  results  on  silicon  [Figure  6-13].  The 
composition  near  the  interface  may  be  slightly  more  Gd-rich, 
which  would  be  consistent  with  a predominance  of  the 
silicon-deficient  orthorhombic  phase,  as  indicated  by  the 
XRD  data  for  the  GaAs  samples. 

Both  Gd  metal  and  GdSi1  7 films  were  deposited  onto  GaAs 
and  interface  stability  was  studied  upon  heating  up  to  500°C 
in  vacuum.  Results  for  rare  earth  metals  and  for  rare  earth 
silicides  were  qualitatively  the  same,  with  greatly  reduced 
reactivity  in  the  case  of  the  silicides.  AES  depth  profiles 
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showed  As  enrichment  at  the  interface,  with  a small  amount 
of  Ga  detected  in  the  silicide  [Figure  6-14].  This  is 
consistent  with  the  observed  reactivity  of  rare  earth  metals 
with  GaAs  ,61"65  where  As-rare  earth  bonding  is  detected  at 
monolayer  coverages  with  the  As  trapped  at  the  interface, 
while  Ga  is  observed  to  diffuse  into  the  metal  film  in  a 
metallic  state. 

There  is  also  a slight  buildup  of  silicon  at  the 
interface,  perhaps  promoted  by  the  loss  of  some  Gd  to  the 
interface  reaction  with  As.  Silicon  once  nucleated  would 
serve  as  a sink  for  excess  silicon  in  the  film,  especially 
if  the  orthorhombic  silicide  is  forming  and  rejecting  excess 
silicon  as  evidenced  by  the  AES  profiles. 

Effects  of  Oxygen 

Oxvaen  Behavior  in  the  Thin  Film  Reaction 

All  rare  earth  metals  have  highly  exothermic  reactions 
with  oxygen,  even  compared  to  the  group  IVA  transition 
metals  [see  Table  6-4].  The  ability  of  gadolinium  to  reduce 
Si02  was  evident  in  our  attempts  to  crystallize  cosputtered 
films  deposited  on  quartz.  Hot  stage  XRD  under  flowing  He 
yielded  no  silicide  peaks,  with  complete  oxidation  of  the 
film  as  low  as  400°C. 

AES  can  provide  depth  information  by  either  sputter 
depth  profiles  or  by  SAM  line-scans  of  craters  from 
stationary  ion-beams. 119  Results  for  partially-reacted  thin 
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films  showed  that  0 was  rejected  from  the  silicide  and 
concentrated  in  the  unreacted  rare  earth  metal  [Figure  6- 
15] ; this  suggests  that  the  silicides  should  be  good 
diffusion  barriers  for  O. 

These  profiles  also  confirm  a mechanism  proposed  to 
explain  the  critical-temperature  reaction  kinetics  of  the 
thin  film  reaction.57  Any  thin  surface  oxide  of  silicon 
would  be  reducible  by  the  rare  earth  metal,  creating  a high 
level  of  oxygen  in  the  rare  earth  metal  that  is  further 
increased  by  the  advancing  silicide  interface.  A thin 
barrier  layer  of  oxide  in  the  rare  earth  metal  is  thus 
created,  hindering  the  diffusion  of  silicon  into  the  metal. 
This  transfer  of  oxygen  to  the  rare  earth  metal  is  confirmed 
in  our  samples. 

The  high  oxygen  level  was  capable  of  stopping  the 
reaction,  approaching  a concentration  sufficient  to  nucleate 
GdO  (AHf=  -17  kCal/mole)  but  not  the  higher  oxide  Gd203  (AHf= 
-434.9  kCal/mole)  , which  is  highly  stable118  [see  Table  6-4]. 
The  formation  of  an  interfacial  layer  of  GdO  at  the  silicide 
interface,  promoted  when  silicide  has  begun  to  form,  is  the 
inhibitor  of  silicide  formation  at  lower  temperatures.  This 
oxide  has  a low  heat  of  formation  compared  to  the  silicide, 
so  its  break-up  is  energetically  favored,  and  proceeds  at 
higher  temperatures. 

The  delay  in  the  onset  of  silicide  formation  has  also 
been  attributed  to  the  presence  of  a carbon  contamination 
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layer.58  In  either  case,  a higher  activation  temperature  is 
required  to  dissolve  or  "punch  through"  the  oxide  layer, 
with  the  silicide  reaction  proceeding  to  completion  within 
minutes. 

Behavior  of  Oxygen  in  Codeposited  Films 

Oxygen  contamination  of  the  codeposited  films  was 
frequently  excessive  due  to  repeated  exposure  of  the  sputter 
target  to  the  atmosphere  in  the  process  of  loading  samples. 
The  effect  of  this  contamination  was  most  readily  observed 
in  sheet  resistance  measurements  [discussed  below],  but  it 
can  be  expected  that  contact  properties  were  affected  as 
well.  Changing  the  apparatus  to  allow  sample  introduction 
while  maintaining  a high  vacuum  in  the  target  chamber  would 
be  highly  desirable. 

AES  analysis  also  showed  that  film  oxidation  in  the 
ambient  produced  a surface  oxide  rich  in  silicon,  while  heat 
treatment  under  a reducing  atmosphere  in  either  RTA  or  a 
tube  furnace  produced  a gadolinium-rich  surface  oxide.  The 
gadolinium  silicide  samples  appear  to  be  protected  against 
extensive  oxidation  similarly  to  ErSi2,  which  is  stable 
against  oxidation  at  about  4 00° C. 75 

Electrical  and  Contact  Properties 

Resistivity 

Sheet  resistance  measurements  using  the  four-point 
probe  technique  were  taken  before  and  after  heat  treatment, 
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with  improved  conductivity  for  the  annealed  films.  Using 
thickness  measurements  from  a quartz  deposition  monitor  as 
well  as  profilometry  data  allowed  calculation  of  the 
specific  resistivity  of  the  silicide.  The  sheet  resistance 
of  the  films  was  about  17  ohms/square  when  150oA  were 
deposited;  this  yields  a resistivity  of  250nO-cm  that  is 
consistent  with  reported  values  of  reacted  thin  films  and 
bulk  silicides9,45,56'71  [see  Tables  2-4,  2-5,  above]. 
Resistivity  was  not  a strong  function  of  Si:Gd  ratio,  but 
increased  to  =500pO-cm  with  oxygen  content  of  =10%. 

Resistivity  in  metal  silicides  has  been  related  to  bond 
lengths  in  the  silicide,5  with  lower  resistivity  observed 
where  metal-silicon  distances  are  shorter  compared  to  metal- 
metal  or  silicon-silicon  distances.  For  orthorhombic 
GdSi1  7,  these  distances  [Table  6-5]  are  consistent  with 
higher  resistivity.  A larger  grain  size  in  the 
recrystallized  films  would  no  doubt  improve  the  measured 
resistivity.  However,  gadolinium  is  among  the  most 
resistive  of  rare  earth  metals  [see  Table  2-5,  above],  and 
other  rare  earths  may  offer  considerably  improved 
conductivity. 

Contact  Properties — Silicon  (100) 

Dot  contacts  (0.7mm  diameter)  for  I-V  measurements  were 
fabricated  by  mechanical  masks  using  moderately-doped 
silicon  (n0=4xl017)  of  both  n-  and  p-type.  The  Schottky 
barrier  height  on  p-Si  was  =0.70  eV  (n  = 1.35),  in  agreement 
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with  the  previously  reported  value.  For  n-Si  the  I-V  showed 
a very  low  barrier,  consistent  with  the  reported  value  of 
0.37  eV6  [Figure  6-16].  Measurements  of  Schottky  barrier 
heights  for  the  rare  earth  silicides  are  difficult  to  obtain 
from  log (I)  vs.  V data  because  of  the  effects  of  high  series 
resistance  and  observed  high  recombination  currents  in  these 
contacts.10,75 

The  difficulty  in  obtaining  barrier  heights  from  I-V 
measurements  and  high  ideality  factors  for  these  silicides 
has  been  attributed  to  the  formation  of  parallel  diodes  with 
a combination  of  high  and  low  barrier  heights.120,121  In  such 
cases,  regions  of  the  contact  having  low  barrier  heights 
will  dominate  the  conduction  even  if  they  comprise  a small 
percent  of  the  total  contact  area.  With  the  increased 
presence  of  large  oriented  grains  in  the  codeposited  films 
such  orientation  effects  may  be  even  more  predominant  for 
our  samples.  Such  nonuniform  contacts  could  also  explain 
the  above-mentioned  high  leakage  currents  observed  in  rare 
earth  silicide  contacts. 

For  more  highly  doped  substrates  of  both  polarities, 
ohmic  behavior  became  predominant  with  little  or  no  heat 
treatment,  suggesting  that  the  films  reacted  with  the 
sputter-cleaned  substrates  due  to  heating  during  deposition. 
The  substrates  were  not  cooled  and  temperature  was  not 
monitored  during  the  deposition. 


113 


Contact  Properties — GaAs  (100) 

Dot  patterns  for  diode  testing  were  prepared  in  the 
same  fashion  on  GaAs  substrates  of  moderate  doping  levels 
(n0=4xl017)  , with  back  contacts  of  annealed  Au-Ge  on  the  n- 
GaAs  and  large-area  silicide  contacts  to  approximate  ohmic 
behavior  for  the  p-type  samples. 

Contacts  to  both  n-  and  p-type  GaAs  were  initially 
rectifying.  Barrier  heights  were  not  obtained,  with  very 
high  ideality  factors  being  found  (n>2)  for  short  linear 
segments  in  the  Log(I)  vs.  V plots  [Figure  6-17].  Because 
recrystallization  of  the  silicide  on  GaAs  showed  even 
greater  orientation  effects  than  on  silicon,  the  comments 
above  concerning  contacts  of  mixed  barrier  heights  is  even 
more  important  for  the  case  of  GaAs. 

Brief  anneals  (600°C/10  sec.)  produced  large  currents 
with  similar  current  response  in  forward  and  reverse  bias. 
In  particular  for  the  case  of  more  highly  doped  substrates 
the  contacts  appear  to  become  near  ohmic  to  both  n-  and  p- 
type,  although  results  were  uneven.  The  decay  of  the 
barrier  properties  and  increased  reverse  current  with  heat 
treatment  are  similar  to  previous  results  with  Pd  on  GaAs.20 
In  this  case  the  same  arguments  for  field  emission122  might 
be  applied,  but  with  indiffusing  Si  as  the  dopant  and 
without  the  incursions  of  reactive  metal  that  give  a 
geometrical  enhancement  of  the  field  emission  for  the  Pd 


contacts. 
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Several  conflicting  theories  of  Schottky  barrier 
formation  continue  to  be  discussed  in  the  literature.66'70 
With  respect  to  this  work,  the  role  of  interface  compound 
formation  may  be  emphasized  for  two  reasons.  Real  contacts 
and  particularly  ohmic  contacts  depend  upon  a certain  degree 
of  reactive  metallurgy,  with  extensive  interfacial  reaction 
products  that  obviate  the  abrupt  interface  model  presented 
in  studies  directed  at  pure  Schottky  behavior.  Secondly,  in 
order  to  eliminate  surface  oxides  and  contaminants  as  much 
as  possible  the  substrates  were  presputtered  before 
depositing  the  silicide.  Sputtering  is  known  to  produced 
altered  surface  composition123  and  broken  bonds  that  will 
promote  interfacial  reaction. 

It  was  shown  above  that  the  annealed  GdSi1  7/GaAs 
interface  was  As-  and  Gd-rich,  with  small  amounts  of  Ga 
dispersed  in  the  film,  consistent  with  previous  studies  of 
rare  earth  metals  on  GaAs. 61-65  These  studies  showed  strong 
arsenide-like  bonding,  including  the  identification  of  an 
interface  phase  comparable  to  bulk  CeAs  in  photoemission 
studies.62  For  the  GaAs  substrates  in  this  study,  free  As 
existing  on  the  sputtered  surface123  would  also  favor  the 
formation  of  an  arsenide. 

The  properties  of  a GdAs  interfacial  compound  is  an 
intriguing  one.  The  rare  earth  monopnictides  (RX,  R = rare 
earth  metal,  X = N,  P,  As,  Sb)  have  attracted  interest 
because  of  their  electrical  and  magnetic  properties.124'125 
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Calculated  band  structures  for  the  gadolinium  monopnictides 
indicate  very  low  effective  masses  for  holes  (0.1mQ  to 
0.3mo)  and  electrons  (0.1mo  in  the  S and  Z directions).126 

The  contact  properties  of  GdAs  on  GaAs  are  not  known, 
but  they  merit  investigation.  The  effects  of  such  a 
compound  at  the  GdSi1  7/GaAs  interface  may  be  involved  in  the 
contact  properties  and  must  be  considered.  GdAs  melts  above 
2270°C127  and  rare  earth  arsenides  are  formed  quite 
exothermically  (HoAs:  AHf=  -72  kCal/mole118)  [Table  6-6]. 

Overall,  the  results  for  GaAs  have  been  inconsistent 
and  unenlightening.  The  effects  of  incorporated  oxygen  in 
the  films  may  also  have  a more  serious  effect  here  than  in 
the  case  for  silicon.  However,  the  possibilities  to  be 
realized  from  a highly  oriented  refractory  film  with  an 
interface  of  limited  reactivity  with  respect  to  the  GaAs 
suggest  that  further  work  is  warranted. 

Problems  in  Patterning  and  Probing 

Electrical  measurements  of  rare  earth  silicide  contacts 
formed  by  codeposition  are  severely  hampered  by  the 
difficulty  in  patterning  the  silicide.  Once  formed,  the 
silicide  is  more  resistant  to  chemical  etches  than  the 
substrate  material.  As-deposited  amorphous  mixtures  with  a 
low  ratio  of  Si/Gd  were  found  to  be  soluble  in  HN03 
solutions,  but  stoichiometric  films  near  the  correct  1.7:1 
ratio  were  not  attacked.  Deposition  from  a sputter  source 
subjects  photoresist  layers  to  polymerization  by  energetic 
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particles  and  heat,  while  the  wide  angular  distribution  of 
flux  from  a large-area  target  does  not  produce  the  clean 
edges  needed  for  effective  liftoff.  For  this  reason,  we 
were  unable  to  produce  suitable  patterns  for  measuring 
contact  resistance. 

The  silicide  is  also  difficult  to  probe  electrically, 
due  to  its  hardness  and  the  presence  of  a surface  oxide. 
Overlayers  of  softer  conductive  metals  helped  with  this 
problem,  but  introduced  more  variables  to  the  guestion  of 
contact  properties. 

Optical  Properties 

The  filling  of  the  shielded  4f  electron  shells  in  the 
lanthanide  sequence  leads  to  a number  of  anomalous 
properties  of  rare  earth  compounds,  including  in  some  cases 
optical  absorption  spectra  that  exhibit  discrete  lines  that 
are  more  typical  of  atomic  or  gaseous  absorption  than  the 
band  absorption  of  solids.2  Thin  layers  of  Yb  metal  have 
been  investigated  as  transparent  top  contacts  for  MIS  solar 
cells.128  It  was  therefore  reasonable  to  investigate  for 
potentially  advantageous  optical  properties  of  the  rare 
earth  silicides,  studies  of  which  have  not  appeared  in  the 
literature  to  date. 

The  optical  reflectance  of  reacted  GdSi1  7 on  Si 
substrates  was  measured  and  a drop  in  reflectance  for 
wavelengths  near  450nm  was  detected  [Figure  6-18].  This 
suggests  that  the  GdSi1  7 has  some  optical  transmission  in 
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the  near  ultra-violet  and  blue  region  of  the  spectrum.  This 
is  further  substantiated  by  the  fact  that  optical 
interference  bands  have  been  observed  for  1000A  GdSi1  7 films 
on  Si  where  the  layer  was  thinned  by  argon  ion  sputtering 
for  AES  depth  profiling  [Figure  6-19].  The  interference 
band  moved  with  the  wavelength  of  monochromatic  light  and 
could  be  seen  for  wavelengths  as  long  as  the  green  region. 
The  transparency  observed  may  be  assisted  by  incorporation 
of  residual  oxygen  in  the  silicide. 

Optical  transmission  measured  on  semi-opaque  400A  and 
80oA  films  indicate  that  any  useful  transparency  is  confined 
to  very  thin  and  hence  very  resistive  layers.  Together  with 
the  reflectance  data,  these  data  are  consistent  with  an 
optical  absorption  length  on  the  order  of  40oA. 

Ellipsometry  measurements  (A  = 6328A)  yielded  real  and 
imaginary  components  of  the  index  of  refraction  similar  to 
NiSi2129 


[Table  6-7] . 
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Table  6-1  X-ray  diffraction  peak  position  (20)  and  d-spacing 
for  5 , OOoA  GdSi2  film  after  various  heat  treatments. 


— 
Sample  treatment 

"1 

20, 

d [ A] 

— 
2 02 

d [ A] 

As  deposited 

34.2 

2.621 

64.9 

1.437 

500°C/20  sec. 

(RTA) 

33.9 

2.644 

64.1 

1.453 

600°C/20  sec. 

(RTA) 

34.1 

2.630 

64.4 

1.447 

700°C/20  sec. 

(RTA) 

34 . 3 

2.614 

65.0 

1.435 

800°C/20  sec. 

(RTA) 

34.4 

2 . 607 

65.5 

1.425 

600° C/40  sec. 

(RTA) 

34.3 

2.614 

64.7 

1.44 

800°C/40  sec. 

(RTA) 

34.4 

2.607 

65.3 

1.429 

375°C/30  min. 

furn. 

34.4 

2 . 607 

64.7 

1.44 

GaAs , 3 7 5° C/ 30 

min. 

34.0 

2.637 

64.1 

1.449 

Table  6-2 
calculated 
(1. 7^x52) , 
references 


Equivalent  lattice  planes  and  d-spacings, 
for  the  three  reported  structures  of  GdSix 
using  lattice  parameters  from  various  cited 
(see  Table  4-6) . 


1 

Crystal  type 

. 

hkl  d [ A] 

hkl  d[A] 

ref. 

Tetragonal 

(112)  2.67 

(220)  1.45 

83 

Orthogonal 

(112)  2.63 

(220)  1.432 

86 

2.626 

1.427 

80 

Hexagonal 

(111)  2.615 

(112)  1.420 

87 

2.606 

1.414 

80 

Table  6-3  Heats  of  formation  for  silicides 


Silicide 

AH./mole 

[kCal] 

ref. 

NiSi2 

-20.7 

4 

WSi, 

-22.5 

4 

TaSi2 

-24.0 

4 

CoSi2 

-24.6 

4 

MoSi2 

-27.9 

4 

TiSi2 

-32.1 

4 

LaSi2 

-43.8 

4 

USi2 

-31.2 
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Table  6-4  Heats  of  formation  for  oxides. 


Oxide 

AH,/mole 
[JcCal  ] 

AHf  per 
oxygen 

ref. 

Sc-,0, 

-411.0 

-137 

4 

w 

-455.5 

-152 

4 

La2°3 

-428.5 

-143 

4 

Gd203 

-434.9 

-145 

118 

A12°3 

-400 

-133 

118 

Zr02 

-261.6 

-131 

4 

Ti02 

-219.0 

-110 

4 

Si02 

-206 

-103 

118 

Ta2°5 

-499.8 

-100 

4 

NbP0, 

-462.7 

-92.5 

4 

wo3 

-200.8 

-66.9 

4 

Cr03 

-136.8 

-45.6 

4 

As205 

-221.05 

-44.2 

118 

Ga203 

-260.3 

-86.7 

118 

Ga20 

-85 

-85 

118 

120 


Table  6-5  Bond  lengths:  (a)  Calculated  [this  work]  for 

orthorhombic  GdSi1  7;  (b)  ThSi2;  (c)  silicon-silicon  bonds, 

related  disilicides. 


(a) 


GdSi2  (orth.) 
Bond  type: 

length  1 

[A] 

length  2 

[A] 

metal-metal 

2.71 

metal-silicon 

2.95 

3.07 

silicon-silicon 

2.22 

2.39 

(b) 


ThSi2 

Bond  type: 

Tetrag. 

i,  [A] 

Hexag. 

1,[A] 

ref. 

metal-metal 

4.135 

4.126 

82 

metal-silicon 

3.161 

3.155 

82 

silicon-silicon 

2.39 

2.31 

85 

(c) 


Silicide 

Tetrag. 

If  [A] 

Hexag. 

1/ [A] 

ref. 

USi2 

2.28 

2.225 

85 

PuSi2 

2.26 

2.24 

85 

YSi2 

2.23 

2.225 

85 

Table  6-6  Heats  of  formation  and  melting  temperatures,  rare 
earth  arsenides. 


AHf/mole 

Arsenide 

[kCal] 

ref. 

HoAs 

-72 

118 

GaAs 

-17 

118 

Arsenide 

Tm 

m 

ref. 

GdAs 

>2270C° 

127 

Table  6-7  Real  and  imaginary  coefficients  of  the  index  of 
refraction  measured  for  various  silicides  (A=6328A) . 


Silicide 

n 

k 

ref. 

GdSi2 

2.85 

2.7 

[PW] 

Nisi 

2.89 

2.87 

129 

MoSi2 

4.16 

1.44 

129 

NiSi2 

3.37 

2.44 

129 

CrSi2 

2.98 

2.13 

129 

WSi2 

4.42 

1.68 

129 

Ir3Si5 

4.50 

1.70 

129 

ReSi 

3.68 

1.53 

129 

CoSi2 

2.13 

1.28 

129 

122 


Figure  6-1  X-ray  diffraction  data:  strong  peaks  identified 
with  the  orthorhombic  GdSi1  , are  indicated. 

(a)  2000A  Gd  deposited  on  Si,  heated  400°C  (fully  reacted) 

(b)  Codeposition,  both  amorphous  as-deposited  and  after 
annealing  at  400°C. 

(c)  Powder  pattern,  bulk-prepared  GdSi2 
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Figure  6-2  X-ray  diffraction  data  for  GdSi2  film  after 
rapid  thermal  annealing  (RTA) : 

(a)  As  deposited  (b)  500°C/20  sec.  (c)  700°C/20  sec. 

(d)  800°C/20  sec.  (e)  800°C/40  sec. 

Thicker  films  show  two  very  dominant  features  [Table  4.3  XRD 
data],  one  near  the  familiar  GdSi2  (112)  peak  at  34°  and 
another  near  20=65°  (d=1.44A)  which  could  be  assigned  to 
GdSi2  (220) . The  dominance  of  these  peaks  is  indicative  of 
two  distinct  strong  preferred  orientations  in  the  films. 
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(c)  orthogonal 
(110)  plane 


(a)  orthogonal 
(112)  plane 


(d)  hexagonal 
(112)  plane 


(b)  hexagonal 
(ill)  plane 


Figure  6-3  Juxtaposition  of  the  hexagonal  unit  cell  and  the 
orthogonal/tetragonal  half-unit  cell,  showing  closely- 
related  planes.  _ 

(a)  orthogonal  (112),  (b)  hexagonal  (111)  planes; 

(c)  orthogonal  (110) , (d)  hexagonal  (112)  planes. 
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Figure  6-4  Expanded  view  of  x-ray  diffraction  data  for 
GdSi2  film  after  rapid  thermal  annealing  (RTA) : 

(a)  As  deposited  (b)  500°C/20  sec.  (c)  600°C/20  sec. 

(d)  700°C/20  sec.  (e)  800°C/20  sec.  (f)  600°C/40  sec. 

(g)  8 00° C/ 40  sec. 

For  the  briefer  20  sec.  anneals,  [(a) -(e)],  the  65°  peak  is 
quenched  dramatically  compared  to  the  as-deposited  film, 
then  grows  slightly  with  temperature,  doubling  from  600°C  to 
800°C  but  remaining  smaller  than  as-deposited  (a) . Longer 
anneals  cause  a pronounced  increase  in  the  peak  intensity 
[(f),  (g) ] . The  marked  shift  to  larger  angle  suggests  a 

phase  change  to  hexagonal  GdSi2. 
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Figure  6-5  X-ray  diffraction  data  for  GdSi2  film  after 
rapid  thermal  annealing  (RTA)  for  different  times: 

(a)  600°C/20  sec.  (b)  600°C/40  sec. 

With  longer  time  at  temperature,  the  65°  peak  increases 
markedly  and  nears  the  intensity  of  the  34°  peak,  indicating 
a time  threshold  for  growth  of  the  grains  contributing  to 
the  65°  peak. 
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Figure  6-6  Expanded  section  of  x-ray  diffraction  data  for 
GdSi2  film  after  rapid  thermal  annealing: 

(a)  As  deposited  (b)  500°C/20  sec. 

(c)  7 00° C/ 20  sec.  (d)  800°C/20  sec. 

The  34°  peak  increases  in  intensity  with  temperature,  and 
appears  at  steadily  increasing  angles  after  an  initial  shift 
to  lower  angle  after  heating  at  500°C.  A distinct  satellite 
peak  at  33.3°  (d=2.69A)  in  the  as-deposited  film  is  totally 
quenched  by  the  500°C  treatment.  At  higher  temperatures  it 
reappears  as  a slight  shoulder  on  the  larger  peak,  becoming 
again  a prominent  satellite  at  800°C. 
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Figure  6-7  X-ray  diffraction  data  for  GdSi2  film  after 
furnace  heat  treatments  (375°C,  30  min.). 

(a)  On  a GaAs  substrate  the  XRD  data  are  dominated  by  a 34° 
peak  assigned  to  orthorhombic  GdSi2(112)  , with  a small  peak 
appearing  at  64°  as  a satellite  on  the  shoulder  of  the 
substrate  peak. 

(b)  On  silicon  the  two  peaks  occur  at  larger  angles  and  at 
nearly  equal  intensity. 

(c)  A spectrum  from  a 800°C/20  sec.  RTA  anneal  is  included 
for  comparison,  showing  the  predominance  of  the  34°  peak  for 
briefer  anneals,  and  the  continued  shift  of  the  65°  peak  to 
larger  angles  at  higher  temperatures. 

(d)  After  1 hr.  at  900°C  a 33.3°  peak,  appearing  as  a 
satellite  in  the  800°C  data,  dominates  a broad  amorphous- 
type  pattern.  No  clear  assignment  could  be  made  for  this  d- 
spacing. 
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Figure  6-8  Planar  TEM  photo. 

(a)  Electron  diffraction  produced  a very  spotty  ring  pattern 
(polycrystalline) . (b)  Bright  field  (BF)  and  (c)  dark  field 
(DF)  images  produced  from  the  indicated  bright  spot  in  (a) . 
The  grains  appearing  in  DF  are  imaged  with  lattice  spacings 
of  ~2.sk. 
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Figure  6-9  Cross-sectional  TEM  of  GdSi2  film  (800C°  RTA 
anneal).  The  camera  axis  is  along  the  substrate  (Oil) 
direction. 

(a)  Electron  diffraction  pattern  shows  spotty  rings 
(polycrystalline) . 

(b)  Bright  field  image 

(c)  Dark  field  (DF)  image  from  spot  C,  lattice  spacing  3.37A 
(hexagonal  GdSi2(110) ) . 

(d)  DF  image  from  spot  D,  lattice  spacing  2.6lA  (hexagonal 
GdSi2(lll) ) . 
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Figure  6-10  Cross-sectional  TEM  of  GdSi2  film  (800°C  RTA 
anneal) . The  camera  axis  is  along  the  substrate  (111) 
direction. 

(a)  Electron  diffraction  pattern  shows  spotty  rings 
(polycrystalline)  with  groupings  of  strong  spots. 

(b)  Dark  field  (DF)  image  from  spot  B,  lattice  spacing 

2.6lA. 

(c)  DF  image  from  spot  C,  lattice  spacing  2.00A;  the  same 
grains  appearing  as  in  (b) . 
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SPOTTER  TIKE 


Figure  6-11  Auger  depth  profiles  of  codeposited  GdSix 
films,  x-2  (silicon-rich),  showing  change  to  1.7:1 
stoichiometry  after  annealing. 

(a)  as  deposited;  (b)  500°C/20  sec., 

(c)  800C°/40  sec.  (RTA) ; (d)  375C°/1  hr.  (furnace) . 
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Figure  6-12  SEM  photo  of  GdSi2  film  after  furnace  anneal  at 
900°C  for  one  hour,  showing  a massive  metallurgical  reaction 
that  produced  silicon  whiskers  growing  out  of  ruptures  in  a 
blistering  surface  oxide. 
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Figure  6-13  Auger  depth  profiles  of  codeposited  GdSix 
films,  x=2  (silicon-rich) , on  GaAs  substrates,  showing 
change  to  1.7:1  stoichiometry  after  anneal. 

(a)  as  deposited  (b)  500°C/20  sec.  RTA 
(c)  700C°/20  sec.  RTA 
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Figure  6-14  Expanded  view  of  Auger  depth  profile  of 
annealed  GdSi2  film  on  GaAs , showing  an  As-rich  interface 
with  a small  amount  of  Ga  detected  in  the  silicide.  Excess 
silicon  is  accumulated  in  the  film  above  the  interface  due 
to  the  loss  of  Gd  in  the  interface  reaction. 
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DISTANCE  X 


Figure  6-15  Auger  line  scan  across  the  edge  of  a fixed  ion- 
beam  sputter  crater  in  a partially-reacted  Gd  film  on 
silicon.  Oxygen  rejected  from  the  reacted  (silicide)  region 
accumulates  in  the  remaining  Gd  metal,  leading  to  the 
formation  of  GdO. 
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Figure  6-16  Current-voltage  measurements  of  GdSi2  contacts 

on  silicon. 

(a)  On  n-type  silicon  - ohmic  as  deposited. 

(b)  On  p-type  silicon  - rectifying  as  deposited,  but 
near-ohmic  after  anneal  of  600°C/10  sec. 

(c)  Log (I)  vs.  V plot  for  contact  to  p-silicon. 

Schottky  barrier  height  measurements  from  these 
data  are  not  truly  meaningful  due  to  high  reverse 
current  and  large  ideality  factor  (n=1.35). 
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Figure  6-17  Current -voltage  data,  GdSi2  contacts  on  GaAs . 

(a)  Contacts  to  p-GaAs  begin  to  show  ohmicity  after  anneal 
of  600°C/10  sec. 

(b)  On  n-GaAs,  rectifying  properties  persist  after  anneal. 
Ohmicity  was  attained  for  more  highly  doped  samples. 

(c)  Log (I)  vs.  V plot  for  contact  to  n-GaAs. 

The  lack  of  a significant  linear  segment  and  high 
reverse  current  make  measurements  of  Schottky 
barrier  height  from  these  data  problematical. 
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Figure  6-18  Optical  reflectance  of  reacted  GdSi, , on  Si 
substrates  showing  a drop  in  reflectance  for  wavelengths 
near  450nm. 

(a)  Reflectance  from  Gd  metal  film. 

(b)  Reflectance  from  GdSi1  7 film. 
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Figure  6-19  Reflectance  measurements  of  small  (50  urn)  spot 

size. 

(a)  A series  of  reflectance  vs.  wavelength  data  taken 
by  stepping  along  the  edge  of  a stationary  ion 
beam  sputter  crater  in  a partially-reacted 
multilayer  sample  (see  Figure  5-5)  consisting  of 
Si/GdSi^/Gd/GdSi,  7/0f-Si. 

(b)  Schematic  diagram  of  the  exposed  layers. 

(c)  Expanded  portion  of  the  optical  data  showing  the 
transition  from  substrate  silicon  to  reflectance 
typical  of  GdSi1  7,  comparable  to  the  data  at  the 
right  from  a codeposited  GdSi1  7 film. 


CHAPTER  SEVEN 
CONCLUSION 

For  several  reasons,  the  rare  earth  silicides  are 
promising  as  semiconductor  contact  materials  for  solar  cell 
applications.  They  have  some  of  the  lowest  Schottky  barrier 
heights  measured  to  date  for  n-silicon.  The  silicides  are 
high-melting,  refractory  materials  with  reasonable 
electrical  conductivity  that  form  reasonably  stable 
interfaces  with  silicon  and  GaAs  substrates.  Ohmic  contacts 
to  n-Si  are  readily  achieved,  and  the  potential  to  form  such 
contacts  with  epitaxial  silicides  is  an  attractive  one. 

Several  modes  of  rare  earth  silicide  deposition  have 
been  investigated,  leading  to  a capability  of  sputter 
depositing  a stoichiometric  mixture  of  Gd  and  Si  with  a low 
level  of  incorporated  impurity,  which  recrystallizes  in  a 
disilicide  structure  with  a moderate  annealing.  The  rough 
surface  of  reacted  silicides  has  been  eliminated  by  the 
codeposition  technique,  and  silicide  films  are  now  readily 
deposited  on  non-silicon  substrates.  The  problems  with 
oxygen  contamination  have  been  shown  to  be  manageable, 
although  improved  consistency  in  film  quality  could  be 
achieved  by  changes  in  the  apparatus  to  eliminate  repeated 
venting  to  atmosphere  during  sample  introduction. 
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The  as— deposited  and  heat-treated  films  show  a strong 
orientation  with  evidence  of  the  hexagonal  phase  being 
formed  on  Si(100),  which  has  not  been  previously  observed. 
Intrinsic  optical  and  electrical  properties  of  these 
codeposited  films  have  been  measured,  as  well  as  contact 
properties  on  both  n-  and  p— type  silicon  consistent  with 
those  reported  in  the  literature. 

In  the  course  of  this  study,  two  novel  methods  of  data 
reduction  were  developed  that  greatly  improve  the 
measurement  of  surface  composition  from  electron 
spectroscopies.  Composition  profiles  derived  from  angle- 
resolved  XPS  data  give  new  insight  into  the  effects  of  ion 
bombardment  on  compounds,  and  discerned  layered  surface 
oxides  on  air-exposed  substrates.  In  addition,  matrix- 
corrected  AES  sensitivity  factors  were  developed  which 
facilitate  monitoring  of  the  composition  for  such  mixed 
films . 

Unfortunately,  the  resistance  of  the  silicides  to 
etching  and  resultant  difficulty  in  patterning  has  prevented 
measurement  of  contact  resistance  even  for  the  case  of  n-Si. 
A full  understanding  of  the  contacts,  including  interfacial 
reaction  with  GaAs,  and  optimization  of  contact  properties 
has  not  been  realized. 

These  silicides  may  yet  be  found  useful  in  improving 
the  performance  of  solar  cells.  Their  potential  as  ohmic 
contacts  to  GaAs  has  not  been  fully  developed.  The  strong 
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orientation  effects  observed  indicate  the  potential  for 
improving  contact  properties  by  increasing  the  area  of  well- 
ordered  interface  associated  with  the  oriented  phase.  Such 
well-ordered  and  oriented  silicide  films  may  have  value  as 
barriers  to  oxidation  or  as  diffusion  barriers  for  other 
conductive  or  distributive  metallizations,  thus  extending 
the  lifetime  of  solar  cell  devices.  Demonstrated  problems 
with  oxidation  of  AlxGa.,_xAs  window  layers130  indicate  that 
consideration  of  rare  earth  silicides  as  protective  barrier 
layers  may  be  justified  independent  of  any  role  in  reducing 
contact  resistance. 
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